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A USER'S GUIDE 
TO THE OA366(3T APL WORKSPACE 



by 

F. Russell Richards 
Naval Postgraduate School 
Monterey, CA 93940 

ABSTRACT 

Instructions are given for the use of the APL Public 
Library Workspace, 2 OA3660, which was developed as an aid to 
interactive exploratory data analysis. The OA3660 workspace 
is accessible to all users of the computer time sharing system 
at the Naval Postgraduate School. The workspace contains 
various data analysis functions, data, and complete internal 
documentation. This report provides a primer on APL, documen- 
tation on each function contained in the workspace, examples 
of the use of each function, and program listings. 
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A USER'S GUIDE 



TO THE OA366jI APL WORKSPACE 



by 

F. Russell Richards 
Naval Postgraduate School 
Monterey, CA 93940 



I . INTRODUCTION 

This report provides documentation for the OA3660 APL 
workspace that was developed during the summer and fall quarters 
of 1977 in conjunction with the offering of the OA3660 course. 
Data Analysis, at the Naval Postgraduate School. The OA3660 
workspace contains APL functions, variables, and data sets that 
are useful for interactive exploratory data analysis. The APL 
functions include selected procedures from the STATlOl public 
library [6], and from Donald R. McNeil's textbook. Interactive 
Data Analysis; A Practical Primer [5] (with some modifications). 
The OA3660 workspace also includes functions written by myself 
and students in my class and various utility functions extracted 
from other public library workspaces. Docxamentation for each 
function is contained in this report and in the OA3660 Public 
Library Workspace. 
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The data sets include data contained in examples and exercises 



in the above mentioned book, data from John Tukey's Exploratory 
Data Analysis [8], data generated from test scores in the 
OA3660 class, and data extracted from a few other textbooks. 

The data are cross referenced with the sources in a variable 
named DATAMAP. 

This report is intended for the user who has some experi- 
ence with APL and who is taking a course in Data Analysis or one 
who is already familiar with the basic techniques of interactive 
exploratory data analysis as described in Tukey [8], McNeil [5], 
and Mosteller and Tukey [7]. This report will not attempt to 
explain the analysis techniques themselves, nor will it discuss 
detailed APL concepts. It will give a brief description of the 
construction of data arrays and transformations of data since 
those are key elements in interactive exploratory data analysis. 
It also includes a brief discussion of basic APL commands such 
as logging into the system, loading workspaces, etc. For further 
details on APL and the NPS time-sharing system, CP/CMS, the 
reader is referred to the Naval Postgraduate School Technical 
Note No. 0141-33, APL(CMS) - An Introduction [1]; the IBM report 
#GH20-0906-l, APL\360-OS and APL\360-DOS User's Manual [2]; 

APL-An Interactive Programming Language [3] by Gilman and Rose; 
and APL Programming and Computer Techniques [4] by Katzan. 
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Each function contained in the OA3660 workspace is 
discussed; the function syntax is given; the function 
parameters (if any) are described; APL listings of the functions 
are given; and at least one example is included to illustrate 
the use of each function . An attempt has been made to assure 
that every function in the workspace is completely debugged. 
However, should any problems be experienced, the user is urged 
to notify me of the problems. Users are also encouraged to 
submit to me interesting data and APL functions that are use- 
ful for interactive exploratory data analysis. 
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APL PRIMER 



II . 

This chapter is written primarily for the user who is 
not familiar with the APL programming language. We attempt to 
provide only that material that the user needs to know about 
APL to use successfully the OA3660 APL Public Library. There- 
fore, we discuss in this chapter login procedures, basic APL 
workspace management, error recovery procedures, creation and 
storage of data, data transformations, function syntax, and 
logout procedures. A user conversant with APL should proceed 
to Chapter III. 

A good understanding of APL would enhance the user's 
facility for working with data arrays and performing interactive 
exploratory data analysis. Therefore, we encourage th user to 
seek out more detailed information on APL. References [1,2, 3, 4] 
are all recommended. 

A. Getting Started; Login 

In order to use the OA3660 APL Public Library, the 
user must have access to a computer terminal that has APL 
capability linked to the Naval Postgraduate School IBM 360/67 
computer. There are many different types of remote terminals 
available for use,* and each has its unique features of oper- 
ation. Therefore, the user should check out the operating 

* 

The IBM 2741 terminals with a special APL typing ball, the 
Intertec terminals, and some of the CRT terminals have APL 
capability. In addition, there are several school owned 
portable terminals with APL capability. 
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instructions for each terminal to augment the general instructions 
given here . 

If the terminal is wired directly to the computer, the 

user need only turn it on to access the NPS CP/CMS Time Sharing 

System. If not, the user must link to the computer via an 

acoustic hookup. First, turn the terminal on to the correct 

settings. Then dial the appropriate telephone number for 

* 

connection to CP/CMS. When the shrill audible tone is received, 
place the telephone receiver into the acoustic coupler which 
is connected to the terminal. If all of the terminal settings 
are correct, and if the time sharing system is in operation, 
the user should receive the message "CP-67 ONLINE". At this 
point, the user must log onto the system by typing (with the 
non-APL character set) : 

LOGIN XXXXPYY 

where XXXX is the user's identification number assigned by 
the computer center, P indicates a private user ID (type G 
if a general user without private disk space) , YY is the terminal 
number, and indicates a carriage return.^ If the identifi- 

cation number is valid, the system will request the user's password. 
★ 

The telephone number and the terminal settings depend on the 
terminal being used. For ASCII terminals (nearly all terminals 
at NPS other than the IBM 2741 terminals) dial either x2611 
(with speed setting at 110) or x3025 (with speed setting at 300) . 
For EBCDIC terminals (IBM 2741), dial x2701. 

On some terminals the user must depress simultaneously the keys 
CONTROL and S in place of a carriage return even though there 
may be a key so labelled. We will use to indicate a 

carriage return for all terminals. 
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The assigned password should be typed followed by a carriage 



return 



PASSWORD (Q) 



If the password is incorrect, the system will reject the login 
and ask the user to start over. If the password is correct, 
the system will request the user's four-digit project number 
to be followed by a four-character cost center code. The 
project number must be assigned by the computer center. The 
cost center code is the user's section identifier or department 
code . 



If the project number is acceptable the system will respond with 
a ready message indicating that the user is in the CMS subsystem 
of the CP/CMS Time Sharing System. Otherwise, the user must 
repeat the entire procedure. 

Once the user has entered CMS he should switch on the 
APL character set* and type 



Switching to the APL character set may consist of flipping a 
switch, changing a type ball or impact print set, or issuing 
a sequence of program instructions depending on the type of 
terminal . 



PPPPCCCC 




APL 
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to enter the APL subsystem. The system will then respond with 
the message: 



A*P*L\N*P*S 

LIBRARY DOCUMENTATION SYSTEM... ) LOAD 1 LIBDOC TYPE DESCRIBE.* 
The user is now in APL, and all of its powerful features 
are available to him. He may use APL in the calculator mode 
somewhat as he would use a hand held calculator; he may define 
his own functions to perform a sequence of operations; or he 
may access the public libraries which consist of commonly used 
preprogrammed functions that perform a variety of useful compu- 
tations . 

In the next section we describe briefly the use of APL 
in the calculator mode to generate, manipulate and store data. 

In later sections we describe the use of APL to transform the 
data, and the use of APL Public Library Workspaces. The user 
is directed to references [1,2, 3, 4] for information about writing 
functions in APL, and for more detailed information about primitive 
operations in APL. 



If any symbol other than the right parenthesis ) appears before 
the word LOAD, the system is not properly translating the 
character set into the required APL characters. If this happens, 
the user should type the symbol that appears (represented here 
by 7 ^) followed by the word OFF as follows: 

T^OFF @ 

This will get the user off the system at which time he should 
consult with the computer center staff as to the special 
requirements for his terminal. Do not proceed in APL until 
the character ")" appears before the word LOAD in the message 
typed when APL is entered. 



9 



B . Generating and Storing Data 

Generation of data arrays is very simple in APL because 

one need not be concerned with format. One must, however, become 

familiar with the APL syntax. First is the assignment symbol 

-t- which plays a role somewhat like the = sign in FORTRAN. 

It means to take the expression on the right and assign it to 

the variable named on the left. Suppose, for example, one wants 

* 

to generate a vector of data named X consisting of the four 
observations; 3.15 12.57 8 6.003. One need only type: 

X ^ 3.15 "12.57 8 6.003 @ 

with the decimal typed if needed and with one or more spaces 
serving as delimiters for the separate observed values. (Real 
numbers may be expressed in scientific format by use of the E 
notation. 2.5E5 means 2.5 x 10^ or 250000.) Extra spaces before 
and after the assignment arrow, , are not needed, but they will 
not hurt anything. Additional values can be added to the vector 
by using the catenate operator (the comma) as follows: 

** 

Just about any name can be used for variables as long as the 
first character is alphabetic and the other characters are 
alphabetic or numeric. The character A is also acceptable. 

§ 

One must take care to strike the characters for the APL 
character set. Several of the APL characters are identical 
in appearance on the terminal keyboards to non-APL characters, 
but different results are obtained. Among the characters 
that appear the same are ) ' (-+xf/:;*. 
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X ^ X, 5.76 “4 1.47 @ 



This statement says to catenate the three indicated values to 
the end of the old data vector X and to call the resulting 
vector X. The entire vector can be viewed by typing X 
followed by the 

X @ 

3.15 "12.57 8 6.003 5.76 "4 1.47 

In the above vector, the negative signs are typed using the 
negative symbol (upper shift 2) , not the minus sign (upper 
shift +) . 

If one is entering many data values into a vector he 
may require several lines of input to do so. Additional values 
could be inserted, as above, using the catenate operator. How- 
ever, a simpler way is to type a comma followed by the quad 
symbol □ (upper shift L) and a carriage return at the end of 
a line if more values are to be continued on succeeding lines. 
The next line will automatically begin with the quad symbol 
prompting the user to continue entering values. This is 
illustrated below; 

DATA ^ 4 6 8 12 6 "2 "4 0 0 1 22 36 29 18, □ 
□ : 4 2 12 9 "3 "17 9 9 2 "1 4 
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If a value is typed incorrectly and discovered before 
the is struck, the value can be corrected by backspacing 

-k 

to the incorrect entry and hitting the line feed key. This 
will erase the incorrect value and everything to the right of 
the value. Then type the correct values. 

Data correction or modification after a line has been 
terminated by a carriage return can be accomplished several ways. 
One way is to determine the index of an incorrect entry and 
assign a new value to that specific element of the data array. 

For example, suppose one wants to change the third element from 
8 to 4.81 leaving all the other elements alone. This can be 
done by typing: 

X[3] ^ 4.81 @ 

Multiple corrections can be made simultaneously as follows: 

X[3, 4, 5] ^ 4.81 5.9 "3.26 @ 

The X vector would then contain: 

3.15 “12.57 4.81 5.9 “3.26 1.47. 

The index operator i (upper shift I) is convenient for 

determining the index of a given value (or indices of a set of 

values) in an array of data. If X is the data array and B 

is a set of values in X, then X i B will generate the set of 
_ 

On some terminals the line feed key is labelled LF. 

On others, the attention or break key must be struck. 
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indices of the values of B in X. If a value appears in more 
than one place in X, only the subscript of its first appearance 
is given. For example, if X ^ 3.5 12.57 4.81 5.9 "3.26 ”4 

the index of the value 5.9 is found as follows: 



X \ 5.9 



4 



Let B ^ 12.57 

given by: 



"3.26. 



The indices of B in X are 



X i B 

2 5 

This operation could be useful for altering selected values in 
an array. To illustrate, assume that we want to change the 
values ~12.57 and ~3.26 in X to 80 and 83. The operation 
below will accomplish this change: 

X[X \ "12.57 "3.26] ^ 80 83 

X 

3.15 80 4.81 5.9 83 "4 1.47 

Additional values can be inserted at any position in the 
vector, and values can be deleted from the vector by using the 
take "t", drop and catenate operations. The take 

operation with syntax 



r t X, 



1.47 
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selects the first r elements from the vector X if r is 
positive and the last |r| elements if r is negative. The 
drop operation with syntax 

r + X 

deletes the first r elements of X if r is positve and the 
last |r| elements if r is negative. If |r| is greater than 
the number of elements in X, r + X will insert zeros to the 
right of the elements of X until |r| elements are obtained 
and r + X will result in the empty vector. The examples below 
demonstrate the take and drop operations. 

X 

3.15 "12.57 4.81 5.9 "3.26 "4 1.47 

3 t X 

3.15 "12.57 4.81 

"2 t X 

"4 1.47 

9 t X 

3.15 "12.57 4.81 5.9 "3.26 "4 1.47 3.15 "12.57 

2 + X 

4.81 5.9 "3.26 "4 1.47 

~4 i X 

3.15 "12.57 4.81 



8 + X 
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The take and drop operators can be used in conjunction 
with the catenate operator to edit data vectors. The examples 
below illustrate some ways this can be done. 

Let 

X^12 34 5 6 78 9 10. 

Create a new vector Y consisting of the first four elements 
of X followed by the elements 11 12 and then the last 
three elements of X; i.e. Y = 1 2 3 4 11 12 8 9 10. 



Y ^ (4 


t X) , 


11 


12, 


( 3 t X) 


0 


Y ^ (4 


t X) , 


11 


12, 


(7 + X) 


@ 



Let W 0 "1 ~2 ~3 ~4 *"5. Create a new vector Y con- 

sisting of the first eight elements of X, followed by the niimber 
25, followed by the middle four elements of W. 

Y ^ (8 t X), 25, (1 1 (5 t W) ) @ 

A few comments about the above operation are in order. APL 
always operates from right to left except when parentheses are 
used to override the standard order of operation. Thus, (5 t W) 
is the first operation executed. This results in the vector 

0 “1 ~2 “3 ~4 . The next operation is (1+0 ~1 ~2 ~3 ~ 4 ) 

which gives ~1 ~2 “3 "4. This is catenated to 25 yielding 

25 "1 ~2 ~3 ~4 which is then catenated to (8 f X) giving 

1 2 3 4 5 6 7 8 25 "1 "2 “3 “4. 
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In data analysis one often wants to compare different 
groups of data or to relate one group to another group. Many 
of the functions in OA 3660 require the data arrays be matrices 
whose columns represent the different groups and whose rows 
represent the various observations for each group. Let us see 
how a data matrix can be created. Suppose for example that one 
wants to create the matrix: 

\ 

1 2 

3 4 

5 6 

\ 7 8 / 

The easiest way to create A is to create a vector consisting 
of the eight elements and then reshape the vector into the 
desired 4x2 array. The reshape operator p can be used to 
reshape a vector into any specified size. The matrix A is 
created by the operation; 

A ^ (4,2)p 12345678 @ 

The general syntax for the reshape operation is A (N,M)pB 
where B is a vector, or scalar, M is the number of columns of 
the array A, and N is the number of rows of A. The first M 
elements of B will be used for row 1, the second M for row 2, 
etc. If there are less than M*N elements of B, the elements 
will be repeated from the beginning as many as are needed. If 
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B contains more than M*N elements, only the first M*N will 
be used. Try a few examples to see what reshape does. 



Remember that the reshape operator fills up matrices 
row by row. Suppose that the user entered the data into the 
vector by columns; e.g.B-^-l 3 5 7 2 4 6 8. How can the 
4x2 matrix A be generated? (4,2) pB will not work since 
that would give 




5 7 

2 4 

\ 6 8 / 



However, if one requested (2,4)pB one would obtain 




What we need now is to transpose this to obtain A. The transpose 

★ 

operator is (upper shift 0 overstruck with upper shift / ) . 

This sequence of operations can be combined as follows: 

A ^ 5i (2,4) pB @ 

The operator p is also useful for determining the 
dimension or shape of an arbitrary variable. When used for 
this purpose it is called the shape operator. The syntax is 




Type shift 0 (oh), backspace, and type upper shift /. 
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Table 2,1 illustrates the results of the shape operator when 
applied to selected variables. The first entry is the scalar 3. 
APL considers all scalars to have no diitiension (not dimension 0). 
Thus, pZ is the empty vector. In the next to last case, the 
variable Z is defined to be lO (iota 0) . (This represents 
a second use of the index operator \ where the operator has 
only a right argument. If K is any nonnegative integer, iK 
is the vector consisting of the first K positive integers. 

For example, i3 would be the vector 1 2 3.) The use of 

the index operator in Table 2,1 is a rather special one, APL 
interprets xO as a vector with 0 elements, hence the result 
0 for pZ. This probably is puzzling to the reader, but it is 
really quite handy. The reader should simply commit to memory 
the facts that scalars have no dimension and the empty vector 
has dimension 0. In the last entry of Table 2,1, Z is defined 
as the 1x3 array consisting of 1 2 3. This array appears 

to be identical to the vector Y -f- 1 2 3, but APL makes a 

distinction between the two. 



VARIABLE 


pZ 


Z ^ 3 


- 


Z ^ 2 3 5 


3 


Z ^ (2,2)p 1234 


2 2 


Z ^ lO 


0 


Z ^ (l,3)p 2 3 5 


1 3 



Table 2.1: THE SHAPE OPERATOR 
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The last APL operators discussed in this section are 
the grade up and grade down functions. These two functions, 

ic 

denoted ^ and t, give the indices of the elements in a 
vector in order of magnitude: ascending order for grade up 

and descending order for grade down. This is very useful for 
data analysis since much of our work requires us to sort the 
data. These operations accomplish this easily. Here are some 
examples : 

A ^ 9 12 6 4 ~2 3 7 "5 11 

856437192 

HA 

291734658 

In the grade up of A the first element, 8, says that the 

eighth value of A should be taken first; the second element, 

5, says that the fifth element of A should be taken next, 
etc., to sort A in increasing order. Similarly for grade down 
to sort in decreasing order. The sorted arrays can be written 
as follows: 

A[M] 

"5 “2 3 4 6 7 9 11 12 

a[Ha] 

12 11 9 7 6 4 3 ~2 ~5 

* 

These symbols are upper shift H and G, respectively, over- 
struck with upper shift M. 
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We have only scratched the surface of APL features 
useful for generating, modifying and storing data. The collection 
of operations discussed here will, nevertheless, enable a user 
to do most of the things that he requires to generate, correct, 
modify and store data. In the next section we discuss arithmetic 
operations in APL so that the user can perform calculations and 
transform data arrays. 



C . Arithmetic Operations 

In this section we discuss some of the basic arithmetic 
operations. We restrict attention to those operations most fre- 
quently used in data analysis . 

The symbols used for the basic operations of adding, sub- 
tracting, multiplying, and dividing are the standard ones used 

* 

in mathematics: ( + , -, x, -r) . However, in APL these operations 

are more powerful than their equivalents in most programming 
languages since the operations can be applied to entire vectors 
or matrices. When used with vectors or arrays, the operations 
are applied componentwise. This requires that the vectors or arrays 
have the same size. An exception is made for the case where one 
of the arguments is a scalar. When this happens APL adds, 

★ 

The subtraction sign is upper shift +. Do not confuse this 
with the negative sign which is upper shift 2. 
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subtracts, multiplies or divides the scalar to each element of 
the vector or array. Of course, the operation can be applied 
to two or more scalars as with any other programming language . 
The examples below illustrate the basic operations: 

A ^ 3 5 7 8 @ 

B ^ ~1 ~2 1 2 @ 

A + B Q 
2 3 8 10 

A - B @ 

4 7 6 6 

A X B @ 

“3 "10 7 16 

A T B @ 

“3 "2.5 7 4 

2 X A @ 

6 10 14 16 

3 + B 

2 14 5 

C ^ (2,2)pA @ 

C 

3 5 

7 8 
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D ^ (2,2)pB (C 



1 2 




1 2 
2 3 



C + B 




8 10 

C V 2 Q) 

1.5 2.5 



3.5 4.0 



3x4+5 




27 



12 V 3 X 2 




2 



The last two examples are included to re-emphasize the 
right-to-left order of operations in APL. Thus, 5 is added to 
4 to give 9 which is then multiplied by 3 to yield 27. Similarly, 
in the last case, 2 is multiplied by 3 to give 6 which is then 
divided into 12 to yield 2. Unlike some other languages there 
is no hierarchy of operations other than right to left. Of 
course, parentheses may be used as in algebra to change the 
priority of operation. 

(3x4) +5 
17 

(12 V 3) X 2 



8 
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In addition to the four basic operations, data analysis 
frequently requires power, log, reciprocal, and exponential 
transformations. As with the four basic operations these can 
be applied to scalars, vectors, or matrices. 

Table 2.2 presents the syntax for these operations. 



Trans formation 


Syntax 


X^ 


X * a 


in X 


@X ^ 


1/X 


vX 


X 




e 


*X 



Table 2.2. SYNTAX FOR DATA TRANSFORMATIONS 



Here we give examples of these transformations : 

X ^ 1 2 3 4 5 @ 

X * 0.5 

1 1.41421 1.73205 2 2.23607 

X * 2 @ 

1 4 9 16 25 

X * "3 @ 

1 0.125 0.037037 0.015625 0.008 



^ The symbol for the logarithm is the upper shift 0 (oh) over- 
struck with the asterisk (upper shift P) . 
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9X 

0 0.693147 1.098612 1.386294 1.609438 




1 0.5 0.333333 0.25 0.2 

X * V “3 @ 

1 0.793701 0.693361 0.629961 0.584804 

(Recall that APL operates from right to left so that the operator 
above raises each element in X to the -1/3 power.) 




2.71828 7.389056 20.085537 54.598150 148.413159 

The logarithm of a number N to an arbitrary base B can be 

determined by typing B @ N; e.g., 

10 9 1000 



3 

Three other useful arithmetic operations for scalars, 
vectors, and matrices are the ceiling f/ floor L, and 
absolute value | . All of these can be used with a single right 
hand argument or with both left and right arguments. 

The examples below illustrate the use of these operators: 



Monadic (one argument) 



3 6 



2 6 



2.75 



T2. 75 
1 "3 

L2.75 
0 “4 

12.75 
6 0.08 



6 0.08 3.6 

(Gives smallest integer argument.) 

6 0.08 "3.6 

(Gives largest integer ^ argument.) 

6 0.08 “3.6 

3.6 (Gives absolute value of each element. 
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3 8 



Dyadic (two arguments) 

3 5 r 2 8 

(Gives the maximum of each component.) 

3 5 L 2 8 

2 5 (Gives the minimum of each component.) 

3 I 7 6 8 ”2 

1021 (Gives the remainder after dividing each element by 3.) 

Clearly, the floor and ceiling operators are useful for 
rounding values. One can select any number of significant digits. 

If one wants to round values to the nearest integer, one should 
type L0.5 + N. If, say, four significant digits are wanted, 
one could type (L.5 + N x 1 e 4) f 1E4 . The examples below 
illustrate the use of L for rounding: 

N ^ 0.0835126 12.51877623 1.33333333 5.25 

L .5 + N 
0 13 1 5 

(L.5 + N X 1e4) t 1E4 
0.0835 12.5188 1.3333 5.2500 

Another operator, the reduction operator /, allows 
the operations discussed above to be applied to all the elements 
of a vector or to the rows or columns of a matrix. The syntax is 
* 1 

All remainders of K|N are expressed as positive integers in 
the set 0, iK-1. If a remainder is negative, K is added to 
it . 
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f/A where f can be any of the arithmetic operations discussed 
above and A is a vector or a matrix. The result of this oper- 
ation when applied to a vector A of size n is the scalar 
A[l] fA[2]f ••• fA[n]. When applied to a matrix M having r 
rows and c columns, the result is the vector B of dimension r 
where B[i] = A[i;l] fA[i;2] f ••• fA[i;c] for i = l,2,...,r 
(column reduction) . The matrix M can be reduced over its 
rows by typing f/[l]M. The examples below illustrate the reductio 
operation : 



15 



84 
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+ /I 3 7 4 

(the sum of all the elements) 
x/1 3 7 4 

(the product of all the elements) 
r/1 3 7 4 

(the largest element) 

M ^ (3,2)p3 1 2 6 5 4 

M 



3 1 

2 6 
5 4 

-/M 

2 41 (differences of elements in rows 1, 2, and 3) 

L/M 

124 (minimum values in rows 1, 2, and 3) 
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+/tl]M 

10 1 (sums of elements in columns 1 and 2) 

r/[l]M 

5 6 (maximtim values in coliamns 1 and 2) 

Finally, we describe the inner product operation and 
the matrix inverse. The inner product, like the reduction 
operation can be applied with any general APL operators. The 
syntax is Af.gB, where A and B are vectors or matrices 
(which must satisfy certain size restrictions) and f and g 
are any general APL operators. For vector arguments, A and B 
must be the same size, say n, and the result is given by 
f/A[l]gB[l] A[2]gB[2] ... A[n]gB[n]. For matrix arguments 

the number of columns of A must be the same as the number of 
rows of B. The result is a matrix of size n x m where 
pA = (n,k) and pB = (k,m). The (i,j)th element of the result 
is f/A [i ;1 ] gB [1; j ] A [i ; 2 ] gB [2 ; j ] ... A [i ;k ] gB [k ; j ] . The 
reader should recognize the operation above as matrix multi- 
plication when f is + and g is x. The examples below 
illustrate the inner product. 

X 1 4 3 2 

y ^ 2 3 1 2 

X + . X y 

21 ((1x2) + (4x3) + (3x1) + (2x2)) 

X X . + Y 

336 ((1 + 2) X (4 + 3) X (3 + 1) x (2 + 2)) 
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A 3 2 p X 6 

A 

1 2 
3 4 

5 6 

B^2 3p35 2 31 1 

B 

3 5 2 

3 11 

A + . X B 



7 


7 


4 


(the matrix product of A and B) 


21 


19 


10 




33 


31 


16 





The domino operator 0 (type □, backspace, is used to 
solve for the matrix inverse of a nonsingular square matrix. 
The inverse of a nonsingular square matrix C is found by 
typing 0C. For example, 

C^2 2pll0 2 
C 

1 1 
0 2 

0C 

1 -0.5 

0 0.5 

C + . X gc 

1 0 
0 1 
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D. Workspace Management 



In the previous sections we have seen how data are input 
and variables are transformed. In this section we describe how 
the user manages his APL workspace so that he can save data for 
use from one session to the next; he can load public library 
workspaces; and he can use functions available in the public 
library workspaces. For information about writing functions the 
user should see references [2,3,4]. 

When a user types APL he is put into a clear APL work- 
space. In this clear workspace data can be created and any of 
the operations described in the earlier chapters can be performed. 
If the user wishes to maintain the data or results for future 
use he must save the workspace. This is done by typing 

) SAVE WSNAME 

where WSNAME is an arbitrary name (first letter alphabetic, 
eleven characters or less) that the user selects for the workspace. 
This private workspace can subsequently be loaded into the user's 
active workspace by typing 

)LOAD WSNAME 

This causes the active workspace to be cleared and a copy of 
the named workspace to be written into the active workspace. 
Alternatively, if the user does not want to clear out the contents 
of the active workspace to bring in the named workspace, he can 
type 
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) COPY WSNAME 



This will simply augment the existing contents of the active 
workspace with the contents of workspace WSNAME. However, if 
an existing variable or function in the active workspace has the 
same name as a function or variable in workspace WSNAME the 
latter will replace the former. 

Any modifications, additions, deletions, corrections, 
etc. that the user makes to a copy of a workspace will affect 
the active workspace, but will not affect the permanent copy of 
the workspace maintained on his private files unless the user 
saves the so modified active workspace. This can be accomplished 
by typing ) SAVE WSNAME (or simply ) SAVE if workspace WSNAME 
were LOADed) . 

The user can determine the name of the active workspace 
at any time by typing )WSID (for Workspace IDentification) . 
Similarly, he can change the name of the active workspace by 
typing )WSID NEWNAME. A list of all the user's private APL 
workspaces is obtained by typing )LIB. Entire workspaces can 
be permanently destroyed by typing ) DROP WSNAME. (Be careful 
with this one I ) 

Within an active workspace, a user can obtain a list of 
all of the functions contained in the workspace by typing )FNS. 
The functions will be listed alphabetically so that all functions 
from those beginning with a certain letter, LETTER, onward can 
be obtained by typing ) FNS LETTER. Similarly, a list of 
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variables in the workspace can be obtained by typing ) VARS or 
) VARS LETTER. Typeout can be terminated at any point by 
hitting the BREAK or the ATTN key. Variables and/or functions 
can be deleted from the active workspace by typing ) ERASE LIST 
where LIST is a single fiinction or variable or a list of functions 
and variables to be deleted. Names in the list should be 
separated by one or more blanks. Two other system commands 
allow the user some control over his APL environment. These 
commands control the number of digits printed out and the width 
of a typed line. They are ) DIGITS N and ) WIDTH N, where N 
is the number of digits to be printed or the desired line width, 
respectively . 

In addition to his own private APL workspaces, the user 
also has access to all of the APL public libraries available 
at NPS . These libraries have numbers between 1 and 999 and are 
intended to hold workspaces of general interest. See [1] for 
a list of the public library numbers. The contents of a pi±)lic 
library can be displayed using the )LIB command followed by 
the library number. For example. 
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)LIB 1 



A- DISK 


R/0 




lOFNS 


17.25 


7/09 


PLOTFORM 


17.25 


7/09 


NEWS 


15.54 


8/27 


WSFNS 


17.25 


7/09 


TEXTEDIT 


17.26 


7/09 


FORMAT 


17.26 


7/09 


CATALOG 


15.17 


7/14 


MAILBOX 


12.50 


7/14 


MULTIPLO 


14.06 


7/18 


FILEFNS 


17.18 


8/14 



The list contains the workspace name and the time and date that 
it was last modified. 

A copy of a public library workspace can be put into the 
user's active workspace by typing )LOAD n WSNAME or 
) COPY n WSNAME where n is the library mamber and WSNAME the 
workspace name. (Recall the differences in LOAD and COPY; LOAD 
will first clear the contents of the active workspace.) Selected 
functions and variables from a workspace (either public or 
private) can be copied by typing 

)COPY n WSNAME OBJECT 
or 

) COPY n WSNAME GROUP 
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whers OBJECT is a singls variabl© or function that ths user wants 
to copy into his active workspace and GROUP is a group of variables 
and/or functions that has previously been defined for V7SNAME.* 

When a public library workspace is brought into the 
user's active workspace, documentation can usually be obtained 
by typing the word DESCRIBE. It is somewhat standard procedure 
to document public library workspaces with a DESCRIBE variable 
giving general information and with each function in the work- 
space being documented by a "HOW" variable. Type the name of 
the function followed by HOW with no intervening spaces . For 
example, documentation for the function STEMLEAF in public 
library 2 OA3660 is obtained by typing STEMLEAFHOW. "HOW" 
variables generally describe the function syntax, parameters, 
input/output requirements, etc. Of course, the user can determine 
if a workspace contains this sort of internal documentation 
by typing )VARS. 

When the user has completed his work session he can logout 
by typing )OFF. If he wants to get out of APL but not logout, 
he should type ) OFF CMS. Then, he should switch back to the 
standard keyboard. 



The group structure for the OA3660 workspace is discussed in 
Chapter III . 
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E . Error Recovery 



This section describes recovery procedures that the user 
can employ if he receives error messages during execution 
of functions contained in a workspace. Usually, error messages 
received when functions in pi±>lic library workspaces are executed 
are the result of improper syntax or problems with the shapes 
of the arguments. As soon as an error is detected, execution 
of the function is suspended, the number of the line containing 
the error is typed, a caret is inserted at the position in the 
line that the error occurred, and an explanatory error message 
is printed out. When the function is suspended the values of 
all variables determined up to the point of suspension can be 
obtained by simply typing the variable name. All of this makes 
error discovery and correction quite simple. From the point of 
suspension the user can branch to any line of the function by 
typing -»■ n, where n is the desired line number. If n is 
omitted a branch is made outside of the function. Many types of 
calculations can be performed while a function is suspended, 
including the execution of other functions. However, the function 
will remain suspended until a branch is made as described above. 
Suspended functions tend to clutter up the user's workspace so 
that suspended functions should not be left pending. A list 
of all suspended functions can be displayed at any time by 
typing )SI (state indicator). If several functions are listed 
as being suspended, the user should type as many branch arrows 
(-»■), one per line, as there are asterisks displayed in the list. 
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since the functions contained in the public library 
workspaces have been tested extensively, most of the errors 
that are encountered result from improper function syntax or the 
use of arguments that are not conformable or of improper size. 
The documentation contained in the next chapter shows the proper 
syntax for each function, and gives the requirements placed on 
the arguments. 
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III. THE OA3660 WORKSPACE 



This chapter provides documentation for the public 
library workspace, 2 OA3660. General workspace documentation 
is contained in the DESCRIBE variable and in the lists of 
functions, variables, and groups. Short writeups are given 
for each function in the workspace via "HOW" variables, and 
examples of the use of each function are provided. Data are 
contained in the workspace to provide easy illustration of the 
functions. The data are described in the DATAMAP variable. 
Finally, the actual APL programs are displayed. 



A. General Documentation 

The OA3660 APL workspace is contained in public library 2. 
Therefore, the user must type ) LOAD 2 OA3660 or ) COPY 2 OA3660 
to create a copy of the OA3660 workspace in his active workspace 
area. The functions and variables contained in OA3660 are 
displayed by typing ) FNS and ) VARS as shown below. 



)Fl\iS 

AiMU AiWVA A2h A2hSR BOXPLOT 

COMPAhh COnDBNSE COi'jT I iJG ENCY 

KS LINE LIT LSLINE 

ONEd ONEZ PARTIAL REGRESS SCAT 
STEMLEAF SUi-LSQ TIICE UTCOND 



CHI SQUARE 
CORRELATION 
EEDPOLISH 
SHOURES SPLIT 



COLE RES 
FILL FE'j 
nSTATS UU'' 
STATISTICS 
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)\/ARS 

AGRS ARDtlO'^l 

CtilCK'^TS 
COlVORRStitiOW 
CRI^^IRS DATAMAP 
DRAPPti fjlPSILOR 
LARDARRAS 
MORTALITY 
RUM RUMSU14H0W 

PRECIPITATION 
RSELECT SCALE SCATHOI SHO'.YRES 

st&mleaehou/ testx testy 

Y ARP b R E AK S RID W I DT H 



ANOVAHOY 
CHI SOU ARE ROY 
CONSUMPTIOi'J 
DEATHS UEP 
GRADESl GRADES'! 
LENGTH LINEHOY 
MSTATSHOY 

ONER HOY 



AERHOY A3RSRR0Y bOKtLO 

COAL CODERESROY CO:PAh 

CONTI NG ENCY HO Y C 0 R Rr.LATIO H 0 -J 
DEPTH ObSChIbE DISCOi' 



HYDROPLANTS 
LSLINEHOY 
NDIVK NOIYY 
0NE3R0Y OPTION 
REG R S S S HOY 
HOY SPLIThOY 

TYICERO’Y USPOP 

^INTERCEPT CR 



I'ySr.CTS HSHOy 
MEDPOLISHROY 
Y G AP N 0 h o r 

PmHTTALHO-J 
RE SIDS Ri ViERS 
STATIS 
VOLCA.^ 



l ,-Gj i 

‘c ;'/■/ 
r I 

L. n I L S 

■1 1 S S 
f ■ / 'T 3 



TICS iOv 
) 



Notice that there is a variable callec! DESCRIBE. This variable 
gives general documentation about the workspace. 



9THE UESCRIdE VARIAbLE PROVIDES GENERAL DOCUhE jT ATI ON 
ftPOR THE ENTIRE YORKSPACE. SISPLY TYPE 'DESCRIbE' TO 
ftO STAIN THE INFORMATION . 



uESCRIbE 

0A3LLQ-DAT A ANALYSIS 

DATE: JANUARY 1978 

PROGRAMMER: P. RUSSELL RICHARDS 

THIS YORKSPACE CONTAINS FUNCTIONS AND VARIABLES USEFUL r ■■ 
EXPLORATORY INTERACTIVE DATA ANALYSIS. FOR INFORMATION Oa T .r. 
USE OF The functions type THE FUNCTION NAME FOLLOYED bY 'HOY'. 
for example, type STEMLEAFHOY for DOCUMENT ATION ON THE FUhCTIO.^ 
STE^iLEAF. SOME FUNCTIONS ARE USED ONLY AS SUBPROGRAMS FOR OTh-,n 
FUNCTIONS AND YILL NOT HAVE A 'HOY' VARIAbLE. FOR :OhE DhThIL~u 
documentation Sba THE NPS TECHNICAL REPORT. A JSE'x'S GUIDE TO 7 ; ■ 
.■J/33660 aPL YORKSPACE. BY E. RUSSELL RICHAbDS On uOUALL . 
YCnEIL'S bOOK. InTERACTIV b DATA ANALYSIS. 

THE YORKSPACE ALSO CONTAINS VARIOUS SMALL DATA SETS tun 
ILLUSTRATION OF ThE FUNCTIONS . TYPE DATAmAP TO ObTAI . 
IDENTIFICATION OF THE DATA ARRAYS. 

THE USER can SAVE SPACE IN HIS ACTIVE YORKSPACE bY COnil > 
ONLY THOSE FUNCTIONS AND VARlAbLbS ■jeEDeU TO PeKFORu / ; 
.bCbSSARY TASKS. FUNCTIONS AND VARIABLES ARn CONVS/IE'TLY 

COLLECTED INTO GROUPS YHICH ARE IDENTIFIED IN THk ' HO i ' 

VARIABLES. FOR EXAMPLE. IF THE USER YANTS TO G‘-FyERATE STn.L 
Or EOXPLOT DISPLAYS HE CAN COPY THE DISPLAY GROUP AS FOLLOW. 

)C0PY 2 (9.43660 DISPLAY 

THE GROUP CONTAINS NOT ONLY THE FUNCTIONS BUT ALSO A! i 
REQUIRED PARAMETERS AND DOC U ME NT AT 1 0 N . 
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There is also a variable named DATAMAP which describes the 
data contained in OA3660 and provides a reference to the data source 
This is shown below. 



fl'/'/ift i/ARIABLt. DATAMAF FROVIDbS TrJ FORMATIU'J AbOUT I'-l". 
f^UATA SSTS COUTAIN&D Jh THF ^ORKSPACt, . IT JIIRS Tun 
f^SOUh'Cb OP SAC a DATA ShT ARD TtiP 37 IP OF Inn DATA A-\La{. 
f\3lyiPLT TIPS 'DATAOAP' FUR DATA DOC U '4S 'JT AT T OR . 



OATAOAP 



TdlS PORtLSPACb CORTAINS SpVSRAL DATA 3‘STS TRAT C An bS USPD 
ILLUoTRATP TtiE FUNCTIOiiS OR TO SER\/S AS !tO'iEPORK PRQELEOS. T 
RbiTiRL COLLECTIUDi OF DATA ARRAYS IS CO>rlAi:!ED IS TRP GhOUP DA.. 
DATA. TtiE DATA CAME FROti THE FOLLOWING SOURCES: 



TO 



{DDRAPElU R.R. ADD H. SMITH, APPLIED uEGRESSIOH 



AS A LYSIS. 



SI LEY ADD SONS. 



^2)ntCSEIL, D.R., INTERACTIVE u AT A ANALYSIS, ADUI SOG-G ESLEY . 
{2,) HILLER, I. AND J.E. FREUND, PhOSAbILITY AnD STATISTICS 
FOR ENGINEERS, PRENTIC E-’d ALL . 

(.R) RICHARDS , F.R., CLASSROOM GRADES. 

(S)TUKEY, JOHn, EXPLORATORY DATA ANALYSIS, ADDISON - V ES LEY . 
THE TAbLE DELOS LISTS EACH DATA ARRAY, THE DIMENSION, TuE SOURCE 
{USING author ADDREV IATIONS) , AND PAGE OR CHAPTER NU'lhERS. 



DATA 


pDATA 


SOU RC E 


AGES 


42 


.MC iP.ll) 


CdICKSTS 


14 6 


MCiP.ZO) 


COAL 


49 


TU 


CONSUMPTION 


5 5 


MC{P. 100) 


CRIMES 


50 4 


MC {P.102) 


DEATHS 


5 4 


94) 


discoveries 


100 


MC{P. 121) 


DRAPE RX 


13 4 


DS{P.1">Q) 


DRAPERY 


1 3 


D5' ( P . 1 7 8 ) 


GRADESl 


15 5 


RI 


GRADES! 


15 5 


RI 


HYDROPLANTS 


34 


TU{CH.O) 


INSECTS 


12 6 


MC{P .12,00) 


LAND A R E A S 


48 


liC{P.O) 


MORTALITY 


5 4 


TU{CH.10) 


PRECIPITATION 


69 


MC{P.O) 


PRESRATING 


114 


MC{P. 120) 


RIVERS 


141 


hC{P. 14) 


TESTX 


10 2 


MF ( P . 2 5 3 ) 


TESTY 


1 0 


MF { P.200) 


US POP 


19 2 


TU{CH.S) 


VOLCANO 


219 


TU 


SARPbREAKS 


9 6 


MC{P .20) 



Finally, notice that there are "HOW" variables for most 
of the functions contained in the workspace. (The only functions with- 
out "HOW" variables are utility functions which are used by other 
functions, but which are transparent to the user.) These variables 
provide documentation on the use of the functions. For example, 
documentation on the function SCAT is obtained by typing 
SCATHOW. 



SCATHOW 

SY.JTAX: SCAT ARRAY 

RARAMETRRS : 



(1) WID~ CONTROLS THE HORIZONTAL SIZE 
{DEFAULT=30 CHARACTERS) . 


OF 


THE 


DISPLAY 


(2) DEP~ CONTROLS THE VERTICAL SIZE 
{.DEFAULT = 1B LINES). 


OF 


The 


DISPLAY 


(3) NDIVX, NDIVY’ NUMBER OF UNITS ON 
RESPECTIVELY {DEF AU LT=H ,H) . 

GROUP: RELATIONS, SMOOTH, COMPARISONS 


X’ 


AND 


Y’AXES, 



DESCRIPTION: SCAT PRODUCES A SCATTER PLOT OP THE DATA CONTAINED 

IN ARRAY. THE ARGUMENT ARRAY CAN BE A VECTOR OR A MATRIX PITH 
MANY AS 9 COLUMNS. IP A VECTOR OF SIZE N, THOSE VALUES ALL 
PLOTTED VS. THE INTEGERS 1 TO N \ IF A MATRIX, THE SECOND, THIRD, 
ETC. COLUMNS ARE PLOTTED VS. COLUMN 1 ON THE SAME AXES. 
DIVISIONS ON THE AXES OF THE PLOT ARE NOT EXPLICITLY PRINTr.0, 
EXCEPT AT THE EXTREMES OF THE PLOT. THE USER CAN CONTROL THE 
RESOLUTION OF THE PLOT BY MODIFYING THE PARAMETERS DID, UbP, 
NDIVX AND NDIVY . PRINTING TIME INCREASES DRAMATIC ALLY NITH 
RESOLUTION -, THEREFORE, OTHER PLOT PROGRAMS SHOULD BE USED IF HIGH 
RESOLUTION IS DESIRED. FOR A SINGLE GROUP OF DATA, THE NUMBER 
{2<N<9) HILL BE PRINTED IF N POINTS LIE CLOSE TOGETHER ON THE 
DISPLAY. FOR MULTIPLE PLOTS ON THE SAME DISPLAY, THE LETTER A 
REPRESENTS GROUP 1, B REPRESENTS GROUP 2, ETC. FOR 2 POINTS 
CLOSE TOGETHER THE LETTER NILE BE PRINTED UITH AN UNDERSCORE. 
THE DISPLAY CANNOT HANDLE 3 OR MORE POINTS CLOSE TOGETHER IP 
THERE ARE MULTIPLE PLOTS. 



The format of SCATHOW is followed for every function. The 
function syntax, the user controlled parameters, the group (or 
groups) containing the functions, the subroutines used by the 
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function, and a brief description of the use of the function 
are displayed. 

The functions and variables are grouped into APL groups 
for ease of handling and to conserve space. The group structure 
is described below. 



)GRPS 

COMPAhlSQiyS 

htILATTOl'IS 



DATA display DOC ESSLUFTALS 

SMOOTH STATS T'^O^YAY 



GkIT t 

i 



)GRP COMPARISONS 
CONDENSE CONDENSEHOW 

SCATdOM LINE LINEHON NDIVX 

COMPARE COMPAREHOH FILL 

ANDtiO'^ LSLINE LSLINSHOW 



NUM MISS UTCOND 

NDIVY NID DEP 

NUMSUM NU^'iSUMHOY 



FMT 
DEPTh 
h SELECT AN 



)GRP DATA 
AGES CHICKNTS 
DISCOVERIES DRAPERY 

INSECTS LAN DARE AS 
RIVERS TESTY TESTY 



COAL consumption CRIMES DEAT>iE 

DRAPERY GRADESl GRADES2 HYDROPLANTS 
MORTALITY PRECIPITATION PRES RATING 

USPOP NARPbREAfiS VOLCANO 



)GRF 

SI Em LE AF 

bCXPLOTHCN 

Ch 



DISPLAY 

STEilLEAFnCN 
LENGTH FILL 



NIDTH SCALE LIT 
nUmSV:! Nl'hSUMHCN 



bCXPLCT 
'■IF F 






)GRP DOC 

ANDHOW ANOVAtiON A2RH0W A3RSRH0N BOYPLOTHON 

CRISUU AREHON CODERESHOW COMPAREHOV CONDENSSHON 

CONTINGENCY HON C 0 RRE LAT I ON HOW DATAMAP DESCRIBE KSiiON 

LINEliON LSLINSHOW MEDPOLISHHOW mSTATHON NUMSUMHOW 

ONeHhOW 0NE2H0N PARTI ALHOW REGRESSHON SCATHON SHOW RESHO^ 

SPLITHOW ST ATISTICSHON STEMLEAFHOW TNICEHON 
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)'jRt' t!3 S kJ 'iTI A LS 



And ANOl/A A3h 


A3kSh bUAf'SO'l" 




b 




COMPARk CON URNS E 


CONTInxjE iCY 


C 0 h’ ti c. L A i 


'ZiJh 


o 


DEPTH EPSILON FILL 


fnt Input 


KS 


Lj b t r . ' . 1 


- 


LSlI .Vl* ms DPULlSii 


I'lISS NDIVa 


j ul ✓ Y 


J J l\ . . t t 


' Cl 


N U i jS U 0 N S H 0 N 2. 3 


0 P '/ 0 V A hG RS S S 


hLbZbS 


bC A b.\ 


J A 


S ii 01 R SS SPL T T S T E M 


Y-'SAF SU ,su 


TvICb-, 


JTr -. jb 


'V 7 


C.N bin TATS PAtxT 


r 




l: rr^i c 





as 



)GRP GFIT 
KSHO^ 



CHISi^UAhS 



C HI S'ciU Ahb', H 0 V 



)GRP kb: LA 2 
Lint. LI.iiF.dCl 
AnUtiC-n kCt.LF.rT 
LI.J'lFkCEP'l 



ICiJF 
LS L 1 1; t: 
SCAT 
b':-i 7 



Lb LI ntlHr'n 
rrAikCG Jib 
rk 



kb'^b.brs 

Ub'.P 



i\ r. u t . 
;; bl\/ .< 



: .1' ' V 



j. i 1 



)GkP SMOOTH 
SPLIT C.-Vi'3 OtVSH 

A3RH0W AZHShliO'^ 
NDTVY ilID DSP 



A3kSh Tires SCAT 
OSEZiiOl OiJSHHOl TlTCSdOl 



SPLTTuO/ 
SCAThOv -.■LI/a 



)GRP FT A 

rcti relation 

rCRhb:LATICNHC>'^ 



k. 

CR Si-iT 

:iFTATFllCl 



tiFTATF PAn-.IAL 
PARTlALtiCl 



FT ATI FT IFF 
FTATIFTlFFuC / 



)GRP TIOIAY 

MSDPOLTSH MSDPOLTSHHOW 

RSSIDS SPSILOn ROHMSFFSCTS 
CODSRSSHOl AND ARDROl 

AnOVAROl OPTION INPUT 



LINS LINSHO-J LSLINS LSLI\L.iO-J 
SiiOlRSS CODERSS RSLLbCT SdOl hto.iO-v 
CONTINGENCY CONTInGl '^CY 101 A'lOl A 

SU :nSO 



41 



The user can conserve space in his active work area by 
selectively loading or copying only the required functions or 
groups. We have attempted to anticipate the type of analyses 
that the user will attempt, and the group structure has been 
selected with the objective to combine functions and parameters 
which naturally go hand-in-hand. A specific group can be 
selected by typing 



)COPY 2 OA3660 GROUPNAME 

where GROUPNAME is any one of the groups listed above. Two of 
the groups, ESSENTIALS and DOC, are especially useful. The 
ESSENTIALS group contains all of the functions and default 
values of all of the required parameters . Documentation and 
data are not included. The DOC group contains all of the 

documentation (DESCRIBE, DATAMAP, and "HOW" variables). 

The user may want to load the entire OA3660 workspace, check 
on some of the documentation, and then erase the group DOC 
to make room for user generated data and other functions or 
variables. The DOC group can be erased by typing ) ERASE DOC. 
Similarly, any group can be erased by typing ) ERASE GROUPNAME. 
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The following sections give brief descriptions (the 



writeups contained in the "HOW" variables) and examples of the 
functions contained in the OA3660 library. Each function descrip- 
tion gives the syntax, a list of the parameters, a list of the 
groups which contain the function, a list of subroutines (other 
functions) used by the function, and a narrative about the use 
of the function. The narrative tells what the function does, 
describes any restrictions on the function arguments, and 
indicates how the function parameters affect the output. The 
order of presentation of the functions follows what appears to 
me to be a natural sequence of exposure to tools of data 
analysis . 

The HOW variables are listed alphabetically in Appendix A, 
and the APL program listings are given in Appendix B. 
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B. Function Descriptions and Examples 



} ii Rp 


UIFPLAY 




FTP:. LEAF 


ftemleafhow 


WIDTH FFALt. 


B 0 X P Ll O '. H 0 'W 


LtNGTti FILL 


nU iFUM HUMFU 


Fn 







C i. /I 1 
' . j C' 



t.i-i Lt A F:i Cvj 

F liFI' AX'. F f Ek L E AF V E FT 0 h 

GhOUFx OlFBLAY 



P A r. A: ic.'.L E E F : 

(1) yiUTti- rCAThOLF TuE PIDTii {FHAkArTEFF Pth 
Tdl' UIFPLAY iUEPAULT^lQ) . 

(2) FCALE- VARIEF ‘TtiE DEPTH (FTEk IdTEkVAL) 
UIFPLAY lE L'HITF OF 1, 2, OR 0.5 TIkEF A PCPEh 
FELEFT Aj IDTEuER FEOk 1 TF 3 iDEFAULT=l) . 

FUbPROGRAM'. LIT 



Lh^E) Ft 

FF Tu; 

Ft 10 , 



DEFFRIPTIOiU FTEPiLEAF GEEERATEF A 
VRFTOR OF CBFEKA ATIFeF . THE FVNFTIC 
DATA UFIEG A FFALIPG RFUTIRE bAFED 
DATA VEFTOR, bY VARYIi^G THE FF ALE 
THE FIZE OF THE FTRM, BUT IT WILL I 



Gib: AidOUHT FPFFIFIED, TRY VARIOUF 



WalFh FHOIFE YIELDF THE BEFT DIFPLAY. 
kAHY FhARAFTERF TO BE PRINTED Ob A 
TRUHFATED AH U THE HUMBER OF TRUHFATED 
AT THE EHD OF THE LIEF, IF THE FTEM 
PO.'/ER OF TEH, THE LEAVER MAY FOHTAIH 
IH ADUITIOH TO THE DIGITF 0 
FhARAFTERF A TO - REPREFENT THE 



FTEM AHD LEAF DIFPLAY 
V AU TOM AT IF ALLY FFALEF 
OH THE RANGE Alib FIZE OF 
PARAMETER THE UFER FAR F 
OT hefeffarily be FFALF 
VALUER LIKE 1,2, AHD 3 

0‘JT AI H F 



IF A LEAF F. 
LInF, the LI: 
FHARAFTERF WILL 
IRTERVAL IF !:■ 
HE FHARAFTERF 
I.v THIF 

NUMbERF 10 TO 19. 



BE 



. U 

vih 
? I 






tlA 

D 



KF 



r/\Fh\ 



A 

' ' ft 

'.L ll !' 

b t/ / 

b Y 
F Eb 

O 

i f I 
/;/• 
V/: :V 
/: 

' «* 3 ‘ 

THE 



FFALE 

1 
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f t/CLCAuC 



00 1 


^ b b & 3 7 9 9 




01 1 


0 0 01 3 5 B 7 9 9 




^ 0 1 


0 0 1 1 0 1 2 4 U u 5 b b b 


7 7 8 3 0 J J 


0 3 i 


•0^ 12244 4.^.^ b 6 b 7 


8 9 0 


Ou 1 


01 1 12 3 3 3 3 3 4 4 4 b 7 


3 6 3b 9 9 


0 3 i 


3 3 1 1 2 ^ 2 5 4 4 *:> 5 6 G G 6 o G 7 7 7 9 9 


0 6 ! 


0 0 1 1 4 4 3 5 0 & 6 G 7 7 7 


8 8 9 


0 7 1 


0 0 0 0 1 1 ^ 2 3 3 4 3 b b 6 


7 8 3 8 9 


0 8 ! 


1 i 2 2 2 J 3 3 '3 b 7 9 




0 9 1 


0 001 2 3344 3 b b 7 7 9 




101 


0112233445339 




1 1 1 


0112334309 




1 / 1 


1124445b 




1 ■> 1 


03478 




1 '* 1 


00 




i 5 i 


067 




1 ^ 1 


2 5 




17 1 


29 





1 5 

1 ‘j 1 0 3 3 7 9 



t i' t. L X L h t Li i H ix 

^ C 1 ^bbbbl A A Ab Juun *' * 

C z I 0 l’ 1 1 2 z 2 4 4 4 ‘i. *j u G 6 7 7 S 8 9 b j /i ;; b ^ F c bU‘^ li 
04 10111^ y^:>bbbbhbl uFFFFiau'j 

9 G I G 0 1 1 4 4 5 S 6 D G G 7 ‘7 7 3 8 9 .4 /i ^ G /> /> r’ F FG :i 1 j. i 

Ob ! 1 ^ 23 bo bG'J 9 A AHbF ly FF F t FG :i ii ■' 

1 C I Z 11 2 3 3 h b ?j b b ^3 ii b b ('' U U F G G ^ 

j -t: 1 i Iz4449G2li^'^. iii 

1 4 I 3 ZGuti 

lo\2bC- 
1 3 ! t) /] b u U ' 



G u A b (J 1 G F G /i br x' G • G u F uF G •• u i :i <*. 
n IbG nFAGI^.u GGauF fILl Abi/AYG IxFFU 
n TilF Fhllb :JUr.bbb^ 



-L 

^ \J IX IL >1 



b'U .Dbb 
L.i' 
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bO X t" LOT n 0 t'J 

b 1 li \L' H A I bwA]^ljO± ^ c. *± 0 1\. 

l-AhAi-iF.^tlh : 

LFFG'Pn- rOtVAhCLF Tut. dORlAOuT uL FIZn CF Tilt. uIFFLh'I 
{URFAULT^bd C H Ah AC? t. uF ) . 

GROUu: DISPLAY 



SUL PRO 
b RFC hi 
OF DAT 
V P P R [\ 
hSTERI 
CbOSSR 
I RTERd 
rhossR 

1 . 5 lit 



GRAi'j I 



t j. 



LL 



PTIOm bOXPLOT uRRt.hATFS A bOA PlOT DISPLAY FOh 
A. A Rt:CT AuGULAh UOX WITU Ed US C 0 REES PO d DI uG TO 
SiUAliTiLES IS PRESEuTEU dlTli TaE i-iEUIAR OAhur., 



r /V . XlIISKEhS ARE URA'RU On EACH 
S aiARKIUG T ti ti LONE ST ARD illbtii.S'S 
UAnTILE distance OF THE UUARTILES, 



SIDE 

data 

LATA 



S iOUTLIEhS) ARE TARiLED NITti ClhCLuS And 
TERQUARTILE DISTANCES GET HEAVY CIRCLES. 



OF THE 
VALUES 
VALUES . 
THOSE 



H VuC, 
LO'/En A 
' NiTit 
LOA .tl 
IT Hi \ 
'UTS I Dp. T 
R 0 h E T n 



.. n 
ND 
hi. 

T u 

A 

r: t 

•1 L 

t-i i. 



2 



UOXFLCT VCLCAiJO 



199 



xxxxxxxx 



★ 



xxxxxxxxxxxxxx OOCO C 

2 2 



9 NE can Squeeze tue loxplot in or stupTch it 

fi CH A N G I N G LEHG T d . 



^ r I rr 

‘b l' .1. 



lY 






2 



[sOXh'LCT VCLCAhC 



199 



xxxx| ★ [xxxxx C0008®® 

I I 32 3 A 



Length*-'? o 

boaflot volcano 

2 

199 



xxxxxxxxxxx] * (xxxxxxxxxxxxxx o OO C O « 



i 



* S « 



46 



i^UMSUM ARhAY 



NUS’iFUihiC^ 

SYiri'AXi MUMSUM VECTOR OR 

PARAi^ItTER I 

.‘iI55n lU'MbFR USEU TO lEDICATR RISFIRG VALUE I.. ..ir uAT/i 
ARRAY iuKFAULT-~22S2'i) . 

GnOUt: CORF ARI SORE 



FU uPkOGRArj i FMT 

uFFCRIPTIORi nUAFUR OPERATEF OR EITHER A VECTOR OR 
IR ALL CAFEF, THE RATA ARE TREATED AF A FIRGLE oA 
PRODUCEF A RUAERICAL FUMPiARY A RICH GIVEF THE FAAPLE 
DELETION OF MIFFInG VALUER) ^ THE EIGhTHF, 

FPREADFy ARD THE MIDPOIRTF IR TABULAR FORM. 
i'FEFL'L FOR TESTIRG THE FYMMETRY OF A OATA FET, ARD T 
THE EFFECTIVE’AEFF OF DATA TRARSFORkATlOuF IR PuODUCIaG 



A.) ARtiAY oUT , 
C(l. RU'iFUE 
FILE {APT.'.i. 

‘l ‘I r. 



THE EXTREhEF, 
THIF FL'MAAaY 



? EVALUAT' 
FY''iMF.TkY . 



RUPiFUi; VOLCARO 



RUkERlCAL FUkkAnY 
FAtiPLE FILE - 219 



; 1 »n r* 1 



kIDFTF 
65. OQ 
66.50 
7 2.50 
100.50 



LOD/8/MIR 

1 

37.00 1 

24.00 i 

2.00 1 

T‘:Tr'“»nnnt'-rr»*^r^n* 



XEDIAR 

65.00 



UPQ/ 8/ MAX 



96.00 1 
121.00 1 
199.00 i 



FPiiE A uF 

i 

59.00 i 

97.00 1 
197.00 1 



n 

fl 

n 



THE FTEHLt.AF, THE BOXPLOT , ARD THE rU MERIC 
FUGGEFT THAT THE VOLCARO DIFT RI cUTIOi^ IF P 
FKEWED. Try a LOG T RARE FORM TO FEE IF TtiA 
THE DIFTRIBUTIOR hORE R EARLY FYMMETRIC, 



AL FV ilk At\Y huL 
OF IT IV ELY 
T JILL aA^kc 



RUkFUi-i ® volcano 



RUkERlCAL FUkElARY 



FAl^iPLE file = 219 



kIDPTF 


LOQi 8/ MI M 


HE Ul A N 


4.17 1 


1 


4.17 


4.09 1 


3.61 1 




3.99 1 


3.18 1 




2.99 1 


0.69 1 








UPQ/8/MAX 

4.56 

4.80 

5.29 



FPhEADF 



0.9 5 ! 

j . 6 2 1 

4.60 1 
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H THERE IS HOW A DECREASING TREND IN THE MIDPTS COLUMN. 
fi THAT SUGGESTS THAT THE LOG IS TOO EXTREME A TRANSFORMATION . 
fl TRI SQUARE ROOT. 



NUMSUM VOLCANO*.^ 



NUMERICAL SUMMARY 
SAMPLE SIZE = 219 

ffn»TWifWWprfn?fn!Tn*»fTrTwrrrfrr«?»^?i"rfirwftwfTn*fnpr’^nrfr»f»»irTrfr?-7fffT^rfr»r^’^nr?f'n»^r*-' — 



MIDPTS 


L0Q/8/MIN 


MEDIAN 


UPQ/S/MAX 


SPREADS 


8.06 1 


1 


8.06 1 


1 




7.94 1 


6.08 1 


1 


9.80 1 


3.72 


7.95 I 


4.90 1 


1 


11.00 1 


6.10 


7.76 1 


1.41 1 


1 


14.11 I 


12.69 



Hnnwi|ifHfTWFfnrTr»w»TWiftfwrT?ffrnT?rfr»Trflrwnrjr»'7»T»7WFwrT'i»^rFr'»»r?t?np'-^T-nT-nr»r»rTsi»Fr-^*^r!r?t'-r»f-»-»>- 



fl THAT LOOKS A LOT BETTER. LET US OBTAIN THE BOXPLOT 
fl AND STEMLEAF PLOT OF THE SQUARE ROOT TRANSFORMED DATA. 



BOXPLOT VOLCANO*.^ 

1.4 

14.1 



O O X X X X X X X X X X X X X X X X X X I * I X X X X X X X X X X X X X X X X X X 000 



LENGTH-^5 0 



BOXPLOT VOLCANO*. 5 

1.4 14.1 



O O xxxxxxxxxxxxx| * |xxxxxxxxxxxxx OO 



fl THE BOXPLOT LOOKS A LOT MORE SYMMETRIC. NOW FOR THE' 
ft STEMLEAF PLOTx 
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01 

02 

0 3 

OU 

0.^ 

0 6 

0 7 

08 

09 

10 

1 1 

1 2 

13 

14 



S'Th.sLi'/.F ^CirAijC* . :> 

u 

2 4 4 U G 
0 0 2 2 2 3 5 

0 0 1 4 4 5 5 G G 7 7 7 9 9 9 
00131223344456677889999 
0 0 0 3 2 2 3 4 4 4 5 6 6 6 6 G G 6 6 3 9 9 9 

0 0 0 0 0 1 1 1 1 2 2 2 3 3 3 4 4 5 5 5 5 5 5 5 S 5 7 7 7 7 8 c 

001 3 1 1 1 3 222 2234444 444 555G7'?773 8oo9 
01113113/334555566677778889 

0 0 0 1 1 3 3 2 2 / 3 4 4 5 5 5 b G 7 8 8 9 

0 0 0 1 3 1 2 2 4 5 G 7 7 8 8 

5 5 5 7 8 

14G899 

G3 
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)GRP rOMPARISOnS 
CORUERSE rORDENSEHOR 

PrATliOW LIRE LINEllGR RblVX 

COMPARE rOMPAREHCW FILL 

ARDHOW LSLIRE LSLIREHOW 



RUM MIFF VICORU 

RDIVY WIU DEP 

RUMFUM RUMFUMhOR 



FM? FCAT 
DFPIh UOAP 
RFFLl-r? Ai.b 



ARDHOW 



FYRTAXi X ARD Y 



GROi'PFi RELATIORF, IWOWAY, COMP ARIFORF 

UEFCKIPIIOR i THE FURCTIOR 'AUD^ IF UFFD TO CREATE A REW 
ARRAY CORFIFTIRG OF X AVGMERTED BY Y AF ADDITIORAL COLUlThF. 
ARD Y CAR BE FCALARF, VF.CTORF, OR MATRICEF BUT CARROT BOTH 
FCALARF, IF AR ARGUkERT IF A FCALAR, A COLUMN IF GENERATED 
ELEMENT OF WHICH IF THE FCALAR. IF THE ARGUMENTF ARE VFCTORF 
MATRICEF, THEY MU FT BE CONFORMABLE . 



DATA 
X 
h E 

FhCH 



fi THE FUNCTION 'AND' IF UFEFUL FOhUATA t.RTRY AND MANIPULA-^ 
fl TIQN. 



A’ •<- 3 4 5 
y-^-9 10 12 
X AND Y 

3 9 

4 10 

5 12 

Z"^3 2 p i 6 
X AND Z 

3 12 

4 3 4 

5 5 6 

5 0 AND Z 

50 1 2 

5 0 3 4 

5 0 5 6 

5 0 AND X 

50 3 

5 0 4 

50 5 



fl EXCEPT WHEN ONE ARGUMENT IF A FCALAR, THE ARGUMENTF OF 
fl ARD kUFT BE CONFORMABLE IN FIZE. 



W'*'^ 6 8 
X AND W 

ARGUMENTF OF 
W AND 

ARGUMENTF OF 



10 






AND 


ARE 


NOT 


AND 


ARE 


l\i 1 



conformable. 
CONFORMABLE . 
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COMPAREtiOW 

SYNTAX i COMPARE MATRIX 
PARAMETERS i 

(1) • DEPTHS VERTICAL HEIGHT OF DISPLAY {DEFAULT =20 
LINES) . 

(2) MISS^ NUMBER USED TO INDICATE MISSING VALUES 
{DEFAULT ="99999) 

(3) NGAP^ NUMBER OF HORIZONTAL SPACES BETWEEN THE 
BOXPLOT DISPLAYS {DEFAULT=3) , 

GROUP t COMPARISONS 
SUBPROGRAM^ FILL 

DESCRIPTION i COMPARE OPERATES ON AN N BY K MATRIX TO PRODUCE K 
VERTICAL BOX PLOTS PLACED NEXT TO EACH OTHER TO ALLOW VISUAL 
COMPARISON OP THE CENTERS, SPREADS, AND OUTLIERS OF THE BATCHES 
{COLUMNS) OF THE MATRIX, THE USER MUST FILL UP THE MATRIX SO 
THAT THERE ARE N OBSERV ATIONS FOR EACH BATCH, THE PARAMETER, 

MISS, SHOULD BE USED TO FILL IN MISSING VALUES, THE USER SHOULD 
ASSURE THAT MISS IS DIFFERENT FROM ALL VALID ENTRIES IN THE DATA 
ARRAY , 



DEPTH 



20 

MISS 

"99999 

NGAP 



3 



fl WE WILL DEMONSTRATE COMPARE WITH THE CHICKWT DATA, 
fl SINCE WE USED 0 TO INDICATE MISSING VALUES IN THAT 
fi DATA, WE MUST FIRST ALTER THE PARAMETER MISS, 



MISS^O 
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COMPARE CHICKWTS 

423 

108 

® 



02 



i 



*I 



T 



O 



1 

X T 




X 



1 




X 

i 



T 



X 



O 



n THE DIAGRAM REVEALS RATHER LARGE DIFFERENCES IN THE 
fi MIDSPREADS OF THE CHICKWT DATA, WE SHOULD TRANSFORM 
fl THE DATA TO TRY TO ACHIEVE HOMOGENIETY IN THE SPREADS, 
fi FIRST LETS GET A NUMERICAL SUMMARY OF THE DIFFERENT COLUMNS, 



CONDENSEHOW 

SYNTAX t CONDENSE MATRIX OR R-^CONDENSE MATRIX 
PARAMETERS \ 

(1) NUMm CONTROLS WHAT STATISTICS ARE INCLUDED IN THE 
SUMMARY, NUMeLl,2,5,T\ (DEFAULT=2) , 

NUM=1 GIVES MEDIANS OF EACH COLUMN OF THE MATRIX 
NUM=2 GIVES MEDIANS AND INTERQUARTILE RANGES 
NUM=5 GIVES MIN ,QUARTILES,MAX , AND SIZE 
NUM=1 GIVES MIN , EIGHTHS , MAX AND SIZE 

(2) MISSr NUMBER USED TO CODE MISSING VALUES 
(DEFAULT=~99999) , 

GROUPx COMPARISONS 
SUBPROGRAMS i UTCOND, FMT 

DESCRIPTION \ CONDENSE GENERATES SUMMARY STATISTICS FOR EACH 
GROUP OF DATA REPRESENTED BY THE COLUMNS OF THE ARGUMENT MATRIX, 
SINCE THE ARGUMENT IS A MATRIX, THE USER MUST FILL UP THE MISSING 
VALUES IN THE MATRIX WHENEVER THE GROUPS (COLUMNS) HAVE DIFFERENT 
NUMBERS OF OBSERVATIONS, THE PARAMETER, MISS, SHOULD BE USED TO 
FILL THE MATRIX SINCE THE FUNCTION WILL RECOGNIZE THOSE VALUES AS 
MISSING DATA AND WILL IGNORE THEM IN ALL CALCULATIONS, THE 

RESULT FROM CONDENSE CAN BE USED AS THE ARGUMENT OF OTHER 
FUNCTIONS SUCH AS SCAT OR LINE IF AND ONLY IF NUM=2, 
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2 



1<1UM 



0 



MISS 



CONDENSE 


CHICKWTS 


151.5 


43 


221 


83.5 


248 


72 


328 


33 


263 


83 


342 


102 



fl THE ABOVE DATA ARE THE MEDIANS AND MIDSPREADS FOR E AC n 
fl OF THE 6 BATCHES (COLUMNS) OF CHICKWT DATA, ZERO IAS 
fl ENTERED INTO THE PARAMETER MISS EARLIER. THE MIOSPREADS 
R SHOW QUITE A DISPERSION , FROM AS SMALL AS 151.5 TO AS LAR3E 
fl AS 342. 



fl WE CAN GET MORE INFORMATION ABOUT THE INDIVIDUAL COLiR’RJS 
fl OF CHICKWTS BY CHANGING NUM AS FOLLOWS’. 



NUM^S 

CONDENSE CHICKWTS 



MIN 


LOQ 


MEDIAN 


UPQ 


MAX 


SIZE 


108.00 


136.00 


151.50 


179.00 


227.00 


10.00 


141.00 


175.00 


221.00 


258.50 


309.00 


12.00 


158.00 


199.00 


248.00 


271.00 


329.00 


14.00 


226.00 


307.50 


328.00 


340.50 


423.00 


12.00 


153.00 


242.00 


263.00 


325.00 


380.00 


11.00 


216.00 


271.50 


342.00 


373.50 


40U . 00 


12.00 



fl IN SEEKING OUT AN APPROPRIATE TRANSFORi^jATION TO EVEN OUT 
fl THE MIDSPREADS , IT IS OS': EN USEFUL TO PLOT THE LOG OF T:iE 
fl MIDSPREADS VS. THE LOG OF THE MEDIANS. 

fl IN MANY CASES WE NEED TO LOOK AT /I SCATTER PLOT OF OUR 
fl DATA, OR AT SCATTER PLOTS OF FUNCTIONS OF OUR DATA. 
fl LET US INTRODUCE A EUNCTION SCAT THAT WILL PROVIDE Ttir 
fl NEEDED SCATTER PLOT. 
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SCATHOy 

SYNTAX I SCAT ARRAY 
PARAMETERS i 



(1) WIDr^ CONTROLS THE HORIZONTAL SIZE 
{DEFAULT -20 CHARACTERS) , 


OF 


THE 


DISPLAY 


(2) dept, controls the vertical size 

(DEFAULT=15 LINES), 


OF 


THE 


DISPLAY 


(3) NDIVX, NDIVYr’ NUMBER OF UNITS ON 
RESPECTIVELY {DEF AULT-H ,H) , 


X^ 


AND 


Y-AXES, 



GROUPi RELATIONS, SMOOTH, COMPARISONS 

DESCRIPTION i SCAT PRODUCES A SCATTER PLOT OF THE DATA CONTAINED 
IN ARRAY, THE ARGUMENT ARRAY CAN BE A VECTOR OR A MATRIX .UTH AS 
MANY AS 9 COLUMNS, IF A VECTOR OF SIZE N, THOSE VALUES ARE 
PLOTTED VS, THE INTEGERS 1 TO N; IF A MATRIX, THE SECOND, TRIED, 
ETC, COLUMNS ARE PLOTTED VS, COLUMN 1 ON THE SAME AXES, 
DIVISIONS ON THE AXES OF THE PLOT ARE NOT EXPLICITLY PRINTED, 
EXCEPT AT THE EXTREMES OF THE PLOT, THE USER CAN CONTROL THE 
RESOLUTION OF THE PLOT BY MODIFYING THE PARAMETERS WID, DEP, 
NDIVX AND NDIVY, PRINTING TIME INCREASES DRAMATICALLY VITH 
RESOLUTION ; THEREFORE, OTHER PLOT PROGRAMS SHOULD BE USED IF HIGH 
RESOLUTION IS DESIRED, FOR A SINGLE GROUP OF DATA, THE NUMBER N 
(2SA/S9) WILL BE PRINTED IF N POINTS LIE CLOSE TOGETHER ON THE 
DISPLAY, FOR MULTIPLE PLOTS ON THE SAME DISPLAY, THE LETTER A 
REPRESENTS GROUP 1, B REPRESENTS GROUP 2, ETC, FOR 2 POINTS 
CLOSE TOGETHER THE LETTER WILL BE PRINTED WITH AN UNDERSCORE, 
THE DISPLAY CANNOT HANDLE 3 OR MORE POINTS CLOSE TOGETHER IF 
THERE ARE MULTIPLE PLOTS, 



30 

15 

4 

4 



WID 

DEP 

NDIVX 

✓ 

NDIVY 



X-^\10 

Yl->-2 3 5 6 8 10 11 13 15 18 
y2-<-20 20 17 24 21 26 31 30 37 33 
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SCAT X AND Y1 
RANGE OF Xi Q 10 
RANGE OF Yi 0 20 



SCAT X AND Y1 AND Y2 
RANGE OF Xi 0 10 

RANGE OF Y\ 0 40 



B B 
B B B 

B 

A A 

AAA 

r>'^^niTTnnnr't'»>Tffrf 



B B 

B B 




A 



A 






r 



ft IF THE ARGUMENT OF SCAT IS A VECTOR, SCAT HILL CREATE 
ft ABSCISSA VALUES TO BE THE POSITIVE INTE GERS 1 TO pARGU.:ENT, 



SCAT Y2 

RANGE OF Xi 0 10 

RANGE OF Yi IS 40 



O 



o 



o 
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I 



n BACK TO THE PROBLEM OF RONHOMOGENEOUS MIDSPREADS lU THE 
fl CHICKWT DATA, LET US LOCK AT A SCATTER PLOT OF THE 
fi LOG MIDSPREADS V, LOG MEDIANS, 

n FIRST WE RESET PARAMETER NUM TO 2, SO IT WILL GIVE 
fi MEDIANS AND MIDSPREADS , 



NUM-^2 



SCAT 9C0NDENSE CHICKWTS 
RANGE OF Xi 5 6 

RANGE 3.4 4.8 



O 



o 



o 



o 



o 



n THE SCATTER PLOT SUGGESTS A RELATIONSHIP BETWEEN MIDSPREADS 
n AND MEDIANS, PERHAPS SOME ROOT OF THE DATA WOULD MAKE 
n THE MIDSPREADS MORE HOMOGENEOUS , TRY SQUARE ROOT, 



CONDENSE CHICKWTS* ,S 



12.30368569 

14.86363273 

15.74801575 

18.11001266 

16.21727474 

18.4912651 



1.717184371 

2.851055496 

2.355341653 

0.9195158937 

2.471407191 

2.852064545 



fi THAT IS BETTER {RATIO OF MAX MIDSPREAD TO MIN MIDSPREAD 
fi IS ABOUT 1.5. MORE WORK NEEDS TO BE DONE THOUGH, 
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)GRP Rh:LATIONS 

LINE LINEHON LSLINE LSLINEHOW 
ANDHOW RSELECT SCAT SCATBOW WID 
^INTERCEPT FMT OR 



REGRESS RECRESSHC‘-J A^u 

PEP NDIVX NDin 



LINEHOW 

SYNTAKi LINE ARRAY OR Z^LINE ARRAY 

PARAMETERi 

RSELECT SELECTS THE OUTPUT MATRIX OF HESILUALS. 
RSELECT=1 GIVES ABSCISSA VS. RESIDUALS {DEFAULTED . 
RSELECT *1 GIVES FITTED VALUES VS. RESIDUALS 



GROUPi RELATIONS 
DESCRIPTION i LINE FITS 



A STRAIGHT LINE 



0 A SET 

BY DIVIDING THE POINTS INTO 3 REGIONS AND USING 
THE X AND Y VALUES IN THE OUTER REGIONS TO DETERM 
THE INTERCEPT IS THE MEDIAN OF THE DIFFERENCES Yt 
ALSO COMPUTES THE RESIDUALS AND GIVES AN N BY 2 M 
EITHER THE X^^VALUES VS. RESIDUALS OR THE FIT 
RESIDUALS DEPENDING ON THE PARAMETER, RSELECT. I 
NOT WANT RESIDUALS TYPED OUT, HE SHOULD USE THE S 

Z^LINE ARRAY 

THE ARGUMENT ARRAY CAN BE EITHER A VECTOR OR A 2n 
IF THE ARGUMENT IS A VEf''^-^ 

CONSTRUCTED TO BE THE FIRt 
IS A Ir^COLUhN MATRIX, THE FIRST COLUMN IS TAKEN T 
X^VALUES AND THE SECOND COLUMN THE SET OF Y’^VALUE 
USED IN CONJUNCTION WITH THE FUNCTION SCAT TO PR 
THE INTERCEPT, AND A PLOT OF THE RESIDUALS. TO D 

SCAT LINE ARRAY 



OF iX,Y) PCI 
THE MEDIANS 
INF THE SLOPE 
SLOPRy-X. L 
ATRIX CONTAIN 
TEu VALUES 
F THE USER D 
YNTAX i 



C 



INF 
I NO 
VS . 
OeS 



TOR OF SIZE N, THE 
N POSITIVE INTEGERS. 



COLUMN MATRIX. 
Xr^VARIABLE IS 
IF THE argument 
BE THE SET OF 
. LINE CA', BE 
^DUCE THE SL':PE, 
THIS, ENTER: 
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16 



X^2 34356767598 10 11 12 14 
OX 



SLOPE: 

2 

3 

4 
3 

5 

6 
7 
6 

7 
5 
9 

8 

10 

11 

12 

14 



ARRAY<-X AND Y 
LINE ARRAY 

1.375 Y^TNTERCEPT: 0.0625 
0.1875 
" 0.1875 
" 1.5625 
0.8125 

1.0625 
" 2.3125 
" 5.6875 
" 0.3125 
" 0.6875 

3.0625 
" 0.4375 

3.9375 

2.1875 

" 1.1875 

1.4375 

0.6875 



RSELECT^O 
LINE ARRAY 

SLOPE: 1 . 3-75 Y ^INTERCEPT: 0.0625 



2.8125 


0.1875 


4.1875 


" 0.1875 


5.5625 


" 1.5625 


4.1875 


0.8125 


6.9375 


1.0625 


8 .3125 


" 2.3125 


9.6875 


" 5.6875 


8.3125 


" 0.3125 


9.6875 


" 0.6875 


6.9375 


3.0625 


12.4375 


" 0.4375 


11.0625 


3.9375 


13.8125 


2 .1875 


15.1875 


" 1.1875 


16.5625 


1.4375 


19.3125 


0.6875 



RSELECT^l 
SCAT LINE ARRAY 

SLOPE: 1.375 Y^INTERCEPT: 0.0625 
RANGE OF X: 2 14 

RANGE OF Y: 4 






o 



o 
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y 

3 4 4 5 8 6 4 8 9 10 12 15 16 14 13 20 

py 

16 

LINh y 

SLOPE: 1.090909091 Y ■’I ETEHCEPT : 0.6818181818 



i 


1 . 22'^2'^2?2'’ 


2 


1. 136363636 


3 


0.04545454545 


4 


“0.04545454545 


5 


1 . 863636364 


6 


“ 1 . 2 2 •;? 2 7 2 2 


r>T 


“4.318181818 


8 


“1.409090909 


9 


“1.5 


10 


“1.590909091 


11 


“O. 6818131818 


1 2 


1.22-7 2 7 27 27 


1 3 


1 .136363636 


14 


1 .954545455 


15 


0.9545454545 


1 6 


1 . 863636364 



16 

ARRAY^X AtW Y 
LI 'YE ARRAY 

SLOPE: 1.090909091 Y^INTERCbPT: 0.6818131313 



1 


1.22"2'^2'’2'’ 


2 


1.136363636 


3 


0. 04545454545 


4 


“O. 04545454545 


5 


1 .863636364 


6 


“l.22'^2'’2'^2'’ 


7 


“4. 318181818 


8 


“1.409090909 


9 


"1.5 


1 0 


1.590909091 


1 1 


“O. 6818181818 


12 


1 . 2272'’2'’2'’ 


13 


1.136363636 


14 


“l. 9 54545455 


15 


0 . 9545454545 


16 


1 . 863636364 



SCAT LINE Y 

SLOPE: 1.090909091 Y -I YTbhCEPT : 0.6818181818 

RANGb OP X: 0 20 

RANGE OF Y: ~6 2 



O 



o 






O 



O 



o 



o 
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RFRLF.CT*-! 

Z^LINF. Y2 

SLCRFx 1.8571U2857 Y ERC F.PT •. 17.428571 43 



LSLIFEHC^Y 

FY,r2A'K X LFLhvF ARRAY OR Z^LFLI'HF ARRAY 

GROUP X RELATIOUF 

DEFORIPTIOH X LFLInE UETF.Ri4lNEF THE LEAPT FQUAREF FOLUTIO'-J, 
Y = Y I^FLCPE^X , THE ARGVkEHT ARRAY 0 AH BE A VtrTOR Oh A l^rcLUL,. 
MATRIX, IF THE ARGUMENT IF A VECTOR, THE X-^VALUEF ARE TAKE^j TO 
BE THE INTEGERF 1 TO N NHERE N^p ARRAY, IE THE ARGUkElT IF /5 
2r-.rCLUMN MATRIX, COLUMN 1 CONFIFTF OF THE X^-VALUEF AND COLUMN 2 
THE Y^VALUEF, REFIDUALF ARE DETERMINED AND FTORED IN A lAThIX 
NhOFE FIRFT COLUMN IF THE FET OF X-^VALUEF AND NnOFE FECOnD COLUMN 
IF THE FET OF REFIDUALF, IF THE UFER NANTF TO FUPPREFF Phl'.TOUT 
OF THE REFIDUALF, HF MU FT AFFIGN THE REFULTF OF LFLINE TO !■. 
VARIABLE-, I,E,, TYPE Z*-LFLINE ARRAY, 



SCAT LFLINE X AND Y2 

FLOPS X 1.933333333 Y r-I nT ERC EPT x 16.8 6 666667 
RANGE OF X\ Q ID 
RAnGE of Y: "4 4 



o o 

o o 






o 



ftIF THE UFER NANTF TO REGREFF A REFPONFF VARIABLE ON 
fi A FET OF CARRIEkF, HE CAN UFE THE PROGRAM REGREFF, 
fi THIF LEAPT FQUAREF MULTIPLE REGREFFION PROGRAM IF 
fi VERY POWERFUL AND FLEXIBLE, 
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REGRESFHCW 

SYRTAXx Z*-Y REGRESS X 
PARAMETER t 

LINTERCEPT^ DETERMINES WHETHER OR NOT AN INTERCEPT TER'} 
IS TO BE INCLUDED, LINTERCEPT=1 GIVES AN INTERCEPT 
TERM, AND RINTERCEPT=0 GIVES NO INTERCEPT, (DEFAULT IS 
1 .) 

GROUP \ RELATIONS 
SUBPROGRAMS i FMT AND SCAT 

DESCRIPTION t REGRESS DOES A MULTIPLE REGRESSION ANALYSIS 
RELATING THE DEPENDENT VARIABLE Y TO A SET OF CARRIERS X, THE 
LEFT ARGUMENT Y IS A VECTOR OF SIZE N, THE RIGHT ARGUMENT X IS 
AN N BY K MATRIX CONSISTING OF N OBSERV ATIONS ON EACH OF K 
VARIABLES OR A VECTOR OF SIZE N IF K=l, OUTPUT CONSISTS OF Ah 
ANOVA TABLE, Ri^SQUARE, STD, ERROR, REGRESSION COEFFICIENTS (THE 
FIRST COEFFICIENT IS THE CONSTANT TERM IF LINTERC EPT=1 .) , 

TvSTATISTICS , V ARIANCEnCOV ARIANCE MATRIX, DURBIN -^W ATSON 

STATISTIC, AND A VECTOR OF PREDICTED Y VALUES AND RESIDUALS, 
THERE IS AN OPTION THAT ALLOWS THE USER TO INPUT A VECTOR OF X 
VALUES AND USE THE REGRESSION EQUATION TO FORECAST Y VALUES, THE 
USER CAN ALSO OBTAIN A SCATTER PLOT OF THE RESIDUALS, WnEN 
EXECUTION TERMINATES , THE PREDICTED Y VALUES AND THE RESIDUALS 
RESIDE IN THE N BY 2 MATRIX Z, 

^INTERCEPT 

1 

pDRAPERX 

13 4 

(\LETS REGRESS A RESPONSE VARIABLE DRAPERY ON THE FOUR 
fi CARRIERS IN THE MATRIX DRAPERX, 
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Z^DRAPERY REGRESS DRAPERX 



ANCVA 



SOURCE 


DF 


SUM SQUARES 


MEAN SQUARE 


ERRATIC 


REGRESSION 


4 


2 ,6619E3 


6. 6697F2 


1 , 1 148F2 


RESIDUAL 


8 


4.7864/?! 


5.9830/?0 




TOTAL 


12 


2 .7158£’3 







R SQUAREi 
STD ERROR i 

COEFFICIENTS 
62.4054 

1 . 

0 . 

0 . 

"O. 

DC YOU WANT 
Y 

VARIANCEnCOVARIANCE 



0.9823756203 

2.446007963 

T STATISTICS 
0.8906 

5511 
5102 
1019 
1441 

A PRINTOUT 



2.0827 

0.7049 

0.135 

“0.2032 

OF THE V ARIANCEr,CCV ARIANCE MATRIX? 



MATRIX X 



4. OOOOf + O 


5.0507£’ + l 


5.0603E+1 5.1660E+ 


1 


4.9597ff 


“5.0507£’ + l 


5. 5468F”! 


5.1266£’“l 5.5425F" 


1 


5.0529? 


“5.0603E+1 


5. 1266£’“l 


5.2387£’“l 5.2570£’“ 


1 


5.1213? 


“5. 166 0£’ + 1 


5. 5425£’“l 


5. 2570/?“! 5.6959f:“ 


1 


5.1688? 


“4.9597E+1 


5. 0529F“1 


5.1213£'“l 5.1688/?“ 


1 


5.0275? 


DURBINmWATSON X 


2.052596933 








DO YOU WANT TO 


FORECAST A 


VALUE FOR Y? 







20 
DC 
Y 

RANGE 

RANGE 



30 15 40 
YOU WANT 



TO SCAT RESIDUALS VS, PREDICTED Y? 



OF 

OF 



Xx 

yi 



70 

~4 



120 

4 






O 
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)GRP SMOOTH 
SPLIT 0NE3 ONEH 

AZRHOW A3RSRH0W 
RDIVY WID DEP 



A2R A3RSR TWICE SCAT 

0NE3H0W ONEHHOW TWICEHOW 



SPLIT hCJ 
SCATHC./ ‘4 bin 



ONEZHOW 

SYNTAX X ONEZ VECTOR 
GROUP X SMOOTH 

DESCRIPTION \ ONEZ SMOOTHS A VECTOR OF DATA 
RUNNING MEDIANS OF 3 WITH TUKEY'S ENDrPOINT 



USING 

RULE. 



OhE 
WHEN S 
VECTOR 
VALUES 



(Jf.y) PAIRS^ THE USER SHOULD USE THE ARGUMENT 
VALUES ORDERED ACCORDING TO THE MAGNITUDE OF THE X 

VECTOR^YLIt.X'l 

IT IS ASSUMED THAT THE X VALUES ARE EQi'ISPACED. 
SMOOTHING BY MEDIANS OF 3 CAN BE ACCOMPLISHED WITH AZR. 
CONSISTS OF THE SMOOTHED SEQUENCE. A PLOT OF THE 
SEQUENCE IS GIVEN BY TYPING x 

SCAT ONEZ VECTOR 



PASS OF 
MOCTnING 
TiiE Y 



REP EAT Lb 
OUTPUT 
SMOOTHED 



X 

3237432 15 4536 



ONEZ X 

233443345453 



SPLITHOW 

SYNTAX X SPLIT AZR VECTOR 
GROUP X SMOOTH 

DESCRIPTION X SPLIT DOES ONE PASS AT DIVIDING MESAS {PAIRS OF 
ADJACENT POINTS WITH A COMMON VALUE WHICH IS A LOCAL MAX OR MIN) 
USING TUKEY'S ENDr<POINT RULE. IT IS USED IN CON ■‘^UNCTION WITH 
AZR. 



SPLIT ONEZ X 

233334345453 



AZRHOW 

SYNTAX X AZR VECTOR 
GROUP X SMOOTH 
SUBPROGRAMx ONEZ 

DESCRIPTION X AZR DOES REPEATED SMOOTHINGS OF RU NNING MEDIANS OF 3 
UNTIL THERE ARE NO CHANGES IN THE SMOOTHED SEQUENCE FROM ONE 
ITERATION TO THE NEXT. SEE ONEZHOW . 
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AZR X 

333443344444 



AZRSRHOW 

SYNTAX \ A3RSR VECTOR 
GROUP \ SMOOTH 

SUBPROGRAMS \ AZR^ 0NE3, SPLIT 

DESCRIPTION X A3RSR DOES REPEATED SMOOTHINGS BY RUNNING MEDIANS 
OF 3 FOLLOWED BY SPLITTING MESAS AND REPEATING UNTIL CONVERGENCE, 
SEE A3R AND SPLIT, 



AZRSR X 

333334444444 



ONEHHOW 

SYNTAX X ONEH VECTOR 
GROUP x' SMOOTH 

DESCRIPTION X ONEH DOES A SINGLE HANNING OF A SEQUENCE OF DATA, 

IT IS USED IN CONJUNCTION WITH THE OTHER SMOOTHING FUNCTIONS /55 A 
FINAL TOUCHmUP TO A SEQUENCE OF SMOOTHED DATA, THE I^rTH RESPONSE 
IS REPLACED BY 0.25xy[I.,ll + O.SxyCJl + 0.25xy[I + l], SEE AZR AND 
AZRSR, 



ONEH X 

3 2.5 3.75 5.25 4.5 3 



9 7 4.25 4.25 6 



TWICEHOW 

SYNTAX X TWICE VECTOR 
GROUP X SMOOTH 

SUBPROGRAMSx ONEZ, SPLIT ^ AZRSR 

DESCRIPTIONx TWICE SMOOTHS A SEQUENCE OF DATA^ THEN SMOOTHS THE 
RESIDUALS AND ADDS THE SMOOTHED RESIDUALS BACK TO THE SMOOTHED 
DATA TO OBTAIN THE FINAL SMOOTH, THE SMOOTHING IS DONE USING 
AZRSR, 
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ICE X 

33333444U4UU 



LEI UE LOOK A'l COME PLCIE 4 II it ■ ^hE I . t. 

R DMA IHA'I CiiODE IHF. EFFECIE CF F^COItllPu. 

:v'ILL USE TtiE FREE RATING DATA. 



WID^Ab 

DEF^2L 

EC AT FREE RATING 
eAnGF of A' t 0 120 

KAnGE CF Yt 0 100 



O O 0 0 

o 2 ° 

o o 

o o 2 
2 o 2 2 

o o 
o 



2® 2^2° ® 2 2°°® 

0 0 o o 2 2 

o2© 3°2®3 000 00 

o 

2 



o 2 
o o 
o o 0^2 
0 2^00 o 






irrO 0f*»p-''^0»-T*-»' 



- - 2 - - 



f^IT IE HARD TO TELL GUC;i FROh TRIE. LETE EKE IF E 
R HELFE ARY. 



EC AT A3REh FREE RATING 
i\ A L! G E CF X i 0 12 0 

uAilGE OF Y i 20 90 



o2 

0 



22333° 

o 

o o 
2 



23 

2 33° 0 

O O 



§2^ 



o 2 o 

o 

22 0023° 

0 0 2 ° 2 
o 2 

o o < 

2 o 



o 2 

o o 



2 

°322 



o 23 3 0 



O O 
2 



z I 
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(^THAT if a lot FIWOTHER. the PFAKF Ai<tD VALLEYF ARE A LOT 
ft riEAREk iJCiV. AF A FINAL TCUONUP, LETF VFE TNICR TO 
ft F'AOOTti THE REFIDUALF AND ADD bAOR TO THE F..OOTHEu FEQVh 
ft AND THEN UFE CNEH AF A FINAL PCLIFH. 



FCAT CNEH TNIOE PREFRATING 
RANGE OF X i 0 120 

RANGE OF Yi 20 90 



O Z 



2 22 3 o 

o o 

O 0 0 
2 2 



2 

2 3 2 2 

233° ° 2 ° o o 23° 

o o 3 o 2 o 2 

2 o 

O O 2 

2 



233 



o 2 

o 

o 2 

o 0 3 2 2 

o o 



o 

o o 
O 2 O 



22 



— T” n n « 
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)GRP TWOWAY 

MEDPOLISH MEDPOLISHHOW 

RESIDE EPSILON NORMEFFECTS 
CODERESHOW AND ANDHOW 

ANOVAHOW OPTION INPUT 



LINE LINEHOW LSLINE LSLIln.RO.-J 
SHOWRES CODERES RSELEr? SHO.-J RKS riOW 
CONTINGENCY CONTINGENCY RC^ A.^rWA 

SU.GSQ 



MEDPOLISHHOW 

SYNTAXi MEDPOLISH MATRIX OR Z^MEDPOLISH MATRIX 

PARAMETERS t 

(1) RESIDSr, CONTROLS OUTPUT OF RESIDUALS. IF 1, 

RESIDUALS ARE PRODUCED IN A TWO-WAY TABLE; IF 2, 

COMPARISON \/ALUES, iRE<> .^CE) iTV AND RESIDUALS Ahi 
PRODUCED IN A TWO COLUMN MATRIX {DEFAULT=l) . 

(2) EPSILON m PROPORTION BY WHICH THE SUM OF THR 
ABSOLUTE VALUES OF THE RESIDUALS MUST BE REDUCED AT 
EACH ITERATION TO CONTINUE POLISHING iDEFAULT^Q .01) . 

(3) NORMEFFECTS’^ DETERMINES WHETHER OR NOT NORMALIZED 
EFFECTS ARE OUTPUT. THE DEFAULT OF 0 SUPPRESSES OUTPUT 
OF NORMALIZED VALUES. NORMEFFECT S=1 CAUSES NOhMALIZED 
VALUES TO BE PRINTED. 

GROUP f TWOWAY 

DESCRIPTION \ MEDPOLISH ITERATIVELY SWEEPS OUT MEDIANS FROm THE 
ROWS AND COLUMNS OF A TWO-WAY TABLE TO YIELD THE MODEL-. 

OBS = MEDIAN f ROW EFFECT -r COLUMN EFFECT RESIDUAL 
IT YIELDS RESIDUALS FOR TESTING THE ADEQUACY OF THE MODEL iSFr. 
SHOWRES AND CODERES). THE FUNCTION CONTINUES POLISHING UNTIL Ttlc. 
STOPPING RULE CONTROLLED BY EPSILON IS ACTIVATED. THE OVERALL 
MEDIAN, THE ROW r, AND COLUMN ’’EFFECTS , AND THE SUM OF ABSOLUT t. 
VALUES OF THE RESIDUALS ARE GIVEN. THE USER CAN SUPPRESS 
RESIDUAL PRINTOUT BY ASSIGNING THE FUNCTION TO A VARIABLE, I.E., 
BY WRITING \ 

Z*- MEDPOLISH MATRIX 



1 

0.01 

0 



RESIDS 

EPSILON 

NORMEFFECTS 
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DEATHS 




11.7 


8.7 


18.1 


11.7 


26.9 


20.3 


41 


30.9 


66 


54.3 



15.4 


8.4 


24.3 


13.6 


37 


19.3 


54.6 


35.1 


71.1 


50 



Z^MEDPOLISH DEATHS 

315.4 

40.3 

38.25 

37.575 

TYPICAL VALUE t 25.6 

ROW EFFECTSt ”15.3625 "lO.l ”l.4 12.6375 32.875 

COLUMN EFFECTSt 2.7 "3.9 12.625 "3.1375 



nTHE RESIDUALS ARE STORED IN Z, THE FIRST VALUES 
p PRINTED OUT ARE THE SUMS OF THE ABSOLUTE VALUES 
R OF THE RESIDUALS, ITERATION CEASES WHEN TWO 
ft .CONSECUTIVE VALUES ARE WITHIN EPSILON, 



Z 

1.2375 2. 

0.1 0 . 

0 0 

0.0625 "3. 

4.825 "O. 



3625 


7.46 


1 


"3.82 




0.17 


4375 


3.73 


275 


0 



25 


1.3 


5 


1.2375 


5 


"1.7625 


75 


0 




"5.3375 



+ /WIZ 

37.2375 



SHOWRESHOW 

SINT AX \ SHOWRES MATRIX OR SHOWRES MEDPOLISH MATRIX 

SUBPROGRAM t AND 
GROUP X TWOWAY 

DESCRIPTIONt SHOWRES AIDS IN THE ANALYSIS OF RESIDUALS FROM 
MEDIAN POLISHING A TWOmWAY TABLE, THE USER SHOULD SET RESIDS=1 
IN MEDPOLISH IF SHOWRES IS TO BE USED, THE OUTPUT OF SHOWRES IS 
A MATRIX CONSISTING OF THE SYMBOLS Oo»xX DEPENDING ON THE SIZES 
OF THE RESIDUALS, LET R REPRESENT THE RESIDUAL, M THE MIDSPREAD 
OF ALL RESIDUALS, L THE LOWER QUARTILE, AND U THE UPPER QUARTILE, 
THE CODES ARE PRINTED AS FOLLOWSx 
Ox R <Ln 1 . SM 



Ox 
O t 
X I 
Xx 



Lml,5M^R^LnM 
L nM <fi <(/ fM 
U+MiR^Ut-l,5M 

R>U+1,5M 
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SHCWkt:? z 

o o Q o 
o o O o 
o o o o 
0 0X0 

X ° ° 0 



CCUKhRStiOW 

syntax I rCDEHF.S UATHIX OR CCDEFRRS EEUEruSn 

SUbPROGRAM AND 
GhOUFi TWONAY 

DESORIPTICN i CODEhES PRCui'CES A DIAGNCRELC ARRAY tOn t -.hLiZi. > 
TRE nESIDUALS FRCk nEuIAn POLlStiING CE A TYO-PAY TAuLt. '■ <r NSr 
SHOULD SET RESlDS^l In hEDPOLISH WREN OCDEReS IS USED. 7nt 
OUTPUT RIATRIX CONSISTS CE TtiE SY^-iBCLS - <> ■^ DEPENDING 0!, E.:r. Sf^- 

CE TRE RESIDUALS. 

VALUE LELOW LOWER UUAhTlLE 
01 VALUE BETWEEN QUArTILES 
+{ value larger TRAN UPPEn qUAETILES 



CODE RES Z 



fi LETS SEE WRAT HAPPENS WiiEN eeSIuS-2. 



RE SIDS*- 2 



Z*-MEDPOLISR DEATHS 

315.4 
4 0.3 
3 8.25 
37.575 

TYPICAL V ALU El 2 5.6 

ROW E ERECT Si ■l5.3 62 5_ "lO.l ”l.4 12.6375 32.875 

CO LUMN EE SECTS i 2.1 3.9 12.625 "3.1375 
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SCAT Z 

RANGE OF X\ ~10 20 

RANGE OF Yt ~8 6 






fl THERE SEEMS TO BE A RELATIONSHIP BETWEEN THE RESIDUALS 
fi AND THE COMPARISON VALUES, (RE* .^CE) iTV . 



A*-LINE Z 

SLOPEt 0.6216761995 Y^INTERCEPTi "0.1348011389 



ANOVAHOW 
SYNTAXi ANOVA 
PARAMETER i 

OPTION m SELECTS THE TYPE OF ANOVA PERFORMED, OPTION=l 
GIVES A ONE-WAY ANOVA AND OPTION=2 GIVES A TWO^NAY 
ANOVA (TREATMENTS WITH BLOCKING), (DEFAULT=1) 

GROUPt TWOWAY 

SUBPROGRAMSx INPUT, SUMSQ, FMT 

DESCRIPTION \ ANOVA DOES A ONEnWAY OR A TWO-WAY ANALYSIS OF 
VARIANCE DEPENDING ON THE VALUE OF THE OPTION PARAMETER, ANOVA 
WILL INTERACT WITH THE USER TO OBTAIN THE REQUIRED INFORMATION 
AND DATA, WHEN ENTERING DATA, SEPARATE THE DATA POINTS WITH AT 
LEAST ONE BLANK AND USE NO OTHER DELIMITERS, OUTPUT CONSISTS OF 
AN ANOVA TABLE, ESTIMATES OF THE OVERALL MEAN, AND ESTIMATES OF 
THE TREATMENT (AND, IF OPTION=2, BLOCK) EFFECTS, 
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Clj r 



CPTICU 



AUCVA 

tJiJTPH NUMBfJR OF FhFATyiFR^P . 



EtiTFR IFCTOR OF lU’kbER OF OFF. FOR TRrATOEnTF 1 TO u 

J! 

3 3 3 3 

y'ETER 3 OhFKRVATIORF FOR TREATEiERT 1 



45 46 61 

EuTER 3 OoFER/ATICFF FOR TREAT’uERT 2 
Li : 

42 44 50 

E'rTER 3 OBFERl/ ATIORF FOR THE AT HE ET 3 

J I 

36 41 48 

Ei-iTt.R 3 ObFERV ATIOuF FOR TREATREiOT 4 

d : 

49 47 54 



FCUROE 


DF 


AnOVA TAbLE 

c C 

V L 


- C 

: i i. 


T RE ATmEnT 


3 


310.92 


3 6.97 


ERROR 


8 


154.00 


19.25 


TOTAL 


11 


264.92 





R-^F-^UAl'iE = 0.419 

OV r. R ALL ! i E A iJ - 
TREATkERT E FFEOTF 



4 6.06 

1.25 "0.75 4.42 3.92 
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fl LET US TAKE THE SAME DATA AND DO A TWO WAT ANCVA, 



OPTION^! 



ANOVA 

ENTER NUMBER OF TREATMENTS , 

□ i 

4 

ENTER tgC-ZOg OF NUMBER OF OBS, FOR TREATMENTS 1 TO u. 

□ i 

3 3 3 3 

ENTER 3 OBSERVATIONS FOR TREATMENT 1 

□ i 

45 46 51 

ENTER 3 OBSERV ATIONS FOR TREATMENT 2 

Di 

42 44 50 

ENTER 3 OBSERV ATIONS FOR TREATMENT 3 

□ i 

36 41 48 

ENTER 3 OBSERV ATIONS FOR TREATMENT 4 

□ « 

40 47 54 



ANOVA TABLE 



SOURCE 



DF 

3 

2 

6 

11 



MS 

20.22 

102.08 

3.64 



F 

5.56 

28.05 



TREATMENT 



60.67 

204.17 

21.83 

286.67 



BLOCKS 

ERROR 

TOTAL 



RwSQUARE = 0 
OVERALL MEAN 



0.924 



45.33 

2.00 0.00 "3.67 1.67 



TREATMENT EFFECTS 
BLOCK EFFECTS *4.58 



0.83 5.42 
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CCNTINGKl\irYHO'.J 
SYinAXi COl<lY'II\.Gt:hrY I-JATRIX 
GhOUFt TV01-/AY 

DESCkIFTIClU rCUTlNGFRCY TAKFF A TWC-^YAY TABLE An j t'r Tn A 
TEFT OF Il'lDEPElWEECE OF kOWF AtiD OOLUkt^F. 7tli'. FuI-FiFAr 
FTATlFTir AbD ITF DEGREEF OF FREEDOM ARE OUTPUT, 



rONTIRGEl\iCY DEATHS 

CHIr^SQUARE VALUE= 2.920832627 DF= 12 



R 

fl 

R 

R 



OOMPARE THE ABOVE VALUE .VITA 
TABLES FOR 12 DF TO DETERlilBE 
STATING THAT THE RO^'=^ 

CAR BE RE-EOTED AT 



-^hOSE ID raI-S.iUAhE 
IF THE HYPOTHESIS 



» AND OOLU'AnS ARE 1 1\> DEFE N DE ;T 
THE SELECTED SIGNIFIC Ai^Ch LEVEL. 
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NSTATS PARTIAL 
PARTI ALtiC'^ 



)GRP STATS 

CORRELATION OR FMT 

CORRELATIONHC^ PiST ATS ROW 



STATISTICS Ch 

STATISTICSliO ! 



ST ATISTICSROW 



SYirTAlx STATISTICS VECTOR 
GROUP t STATS 
SUBPROGRAM i FMT 

description : STATISTICS DETERMINES THE MEAN, VARIANCE 

DEVIATION, COEFFICIENT OF VARIATION, LOWER AND UPPER 



S X A i’i D A) u 
QUART ILFS, 



MEDIAN, Thl'-'iEAN, MIDMEAN, MIDRANGE , RANGE, MEAN 

DEVIATION, Interquartile range, skewness, and kurtosis 

DATA IN VECTOR, 



ADSCLUTE 
FCh THE 



STATISTICS VOLCANO 
MEAlU 70.24657534 
VARIANCEx 1850.562775 
STD. DEV.i 43.01816796 
COEFF. OF VARIATIONi 0.6123881165 
LOWER QUARTILEi 37 
UPPER QUARTILEi 96 
MEDIANi 65 
T RIME AN i 6 5.75 
MIDMEAN i 64.47747748 
RANGE i 197 
MiDRANGEi 100.5 

MEAN ABSOLUTE DEVIATION i 33.68493151 

Interquartile range i 59 

COEFF. OF SKEWNESSt 0.8325767059 
COEFF. OF KURTOSIS i 0.4281831164 



correlationhow 

SYNTAXi R^CORRELATION W 
GROUPi STATS 
SUBPROGRAM: FMT 

DESCRIPTION t CORRELATION DETERMINES THE SIMPLE PEARSON PRODUCT 
MOMENT CORRELATIONS BETWEEN EACH PAIR OF VARIABLES REPRESENTED BY 
THE C COLUMNS OF W. THE OUTPUT IS A C BY C CORRELATION MATRIX. 
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f ) IN THE rCHhELATlOlJ OUTPUT , THE TCP HOV AHb T : ii -. blHTT 
p COLUkN ARE LAHELE ECh TiiE OC nRE LkT TO N E . FOR >‘' AA .. tLE , 
p THE rCRRELATlCN BETWEEN THE FFrCNu Ai.D ECURTH 
p OF GRAOESl IE 0 , 5 ^, 



pGRADEEl 

15 5 



ESTATE HOW 
EYNTAX : MET ATE U 

GROUP i ETATE 



EUBPRQGRA'Ai i i^iT 



UEECHIFTION \ AETATE UETERAINEE CC RRE L AT 10 E , 
BE \/ IATICNE , LOWER QUARTILEE , AEuIAHE , AnD UPPER 
COLU - inE of TriE ARGUAEnT i-iATRIX W . THE CORRELAT 
IN AN N BY N kA ' TRIX , WHERE N IE THE NUAbER CE 0 



■lE A E , 
N.U A RT I L r. 
IChE ARt 
OLUA.E Ot 



E'T A u L ' kr . 
E 1 0 :x T i 
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/ . 
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PARTIALHOW 
SYNTAX i PARTIAL 
GROUP'. STATS 

SUBPROGRAMS'. OORRELATION , FMT 

DESCRIPTION i PARTIAL DETERMINES PARTIAL CORRELATIONS OF At.Y 
SPECIFIED ORDER FOR THE COLUMNS OF THE MATRIX N. THE USER GILL uE 
ASKED FOR THE COLUMNS TO BE PARTI ALLED OUT. THE USER'S RESPO.^SeS 
MUST BE INTEGERS BETGEEN 1 AND C, GHERE C IS THE NUMBER OF 
COLUMNS OF W. 



PARTIAL GRADESl 

ENTER COLUMuS OF MATRIX TO BE PART I ALLED OUT. 
J: 



1 
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2.0 0 
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fi BE CAREFUL WITH INTERPRETING THE LABELS IN THE OUTPUT 
ft OF PARTIAL. ABOI/E, LABEL 1 DOES NOT MEAN COLUMN ONE, 
ft BUT THE FIRST COLUMN REMAINING AFTER THE SELECTED VARIAbL 
ft HAVE BEEN PARTI ALED OUT. SIMILARLY, 2 REFERS TO TtiE 
ft SECOND REMAINING COLUMN. 



PARTIAL GRADESl 

ENTER COLUMNS OF MATRIX TO BE PARTIAL LED OUT. 
Ut 



2 3 
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ft THE ABOVE EXAMPLE SHONS THAT MORE THAN ONE VARIABLE 
ft CAN BE PARTI ALED OUT. LABEL 2 ABOVE REFERS TO THE 
ft SECOND COLUMN REMAINING AFTER THE SELECTED VARIABLES 
ft HAVE BEEN PART I ALED OUT, IN THIS EXAMPLE 2 REFERS TO 
ft COLUMN k. 
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jGRP GFIT 

KB KBtiCW CdIBQUARK rHIBQL'AhRdCP 



CHIB-QU ARFJWG 
SYliTAXx T CHIBQUARF X 
GROl'Pi GFIT 

DF.BChlPTIOi^ i CHIBQl'ARE CCAPARFF A T RFC RET IC Ah UIBritT" 
PROBABILITY MABB FURCTICN T HAVIEG N PCBBIhLE VALU^B .-JlTH h, 
EMPIRICAL FhEQUERCY FUNCTION X OVLR THE BALE N VALUEB TO TEBT :iC.:' 
l^ELL THE PRObABILITIEB T FIT THE DATA X, THE Cnl-BKi' ARE 
GOODNEBBf^OFnFIT BTATIBTIC IB PRINTED OUT NITH THE NUhBEE, Ci 
DEG REE B OF FREEDOM. 



fl THE THEORETICAL PROBABILITY MABB FUNCTION T A^D THE 
fl EMPIRICAL FREQUENCY FUNCTION X MUBT BE OF THE BAmF 
fl BIZE\ I.E., pf=pA. THE ACTUAL OBBERVED VALUEB ARE 
fi IMPLICIT IN THE PO BIT ION OF THE VECTORB. 



T^.l .2 .3 .2 .1 
X^ 5 20 38 32 12 



THE VALUEB IN T REPREBENT PR 
X REPREBENT OBBERVED FREQUENCIEB . 

ARE NOT BTATED EXPLICITLY , BUT IT 
THE POBITIONB OF X AND T REFER TO THE 
EXAMPLE, IN THE CABE ABOVE, THE FIRBT 



THE VALUEB I.. 
THE ACTUAL OBB'^RVATIO 
IB UNDERBTOOu T.-imT 

BADE VALUEB. FOR 
POBITION MAY REFEh 



rp r 



HE VALUE 



THE BECOND 0, ETC. 



T CHIBQUARF X 

CHI-BQUARF. VALUE = 9.62087227U DF - 4 



fl TO DETERMINE IF ONE BHOULD hF -EOT THE HYPOTHEBIB T n AT 
fl THE DATA HERE DRA.VN FROM A POPULATION IITH PROBABILITY 
fi FUNCTION T, COMPARE THE. VALUE ABOVE UTH THOBF 

fl IN CHI^BgUARE TABLEB pith 4 DE . 
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SYi\ITAXt 

GRCUPt 



c nc 
GFIT 



S G 



DESCRIPTICR i KF DOES A RO Ll40G0RQ\/ ^FklRNCV GOODREFF 



OF TWO OIFCRF.TK 
THE ARGUEiERTF F ARQ 



OF THE EQUALITY 
DIFTRIbUTIORF , 

OBFEkV ATIORF, TREY REED ROT bt FCRTED 
FIZE, OUTPUT OF KF ARE THE VALUEFi 

hAX F{X)vG{X) ARE MIN F{X)-!G{X) 



^ t IT 
EMPIRICAL PHObAb 
G ARE EACH VECTOR 
NOR MU FT THEY bE THE 



TEFT 
iLITY 
F OF 
FAME 



F*-3 259101635719201263889 7 77 7 2 137 
G^H 2 10 20 18 3 6 7 U 2 1 10 8 6 3 31 7 



F KF G 

max FiX)fGU) =0.1274509804 
Mik FiX)r:G(X) ='0.3137254902 



ft 

ft 

ft 

ft 



TO UETERMIR'E IF THE TWO FAMPLEF OF DATA HFRF 
THE FAME POPULATIOR {.OR POPULATIOUF WITH THE 
COMPARE THE MAX ARD MIR VALUEF WITH THE DELTA 
IN TAbLEF FOR THE TWO F AMPLE KF TEFT, 



TAKER FROM 
FAME CDF'F) 
FTATIFTICF 
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APPENDIX A 
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AN DHOW 

SYNTAX: X AND Y 

GROUPS: RELATIONS, TWOWAY, COMPARISONS 

DESCRIPTION: THE FUNCTION 'AND' IS USED TO CREATE A NEW DATA 

ARRAY CONSISTING OF X AUGMENTED BY Y AS ADDITIONAL COLUMNS. X 
AND Y CAN BE SCALARS, VECTORS, OR MATRICES BUT CANNOT BOTH BE 
SCALARS. IF AN ARGUMENT IS A SCALAR, A COLUMN IS GENERATED EACH 
ELEMENT OF WHICH IS THE SCALAR. IF THE ARGUMENTS ARE VECTORS OR 
MATRICES, THEY MUST BE CONFORMABLE. 



ANOVAHOW 
SYNTAX: ANOVA 

PARAMETER: 

OPTION- SELECTS THE TYPE OF ANOVA PERFORMED. 0PTI0N=1 
GIVES A ONE-WAY ANOVA AND 0PTI0N=2 GIVES A TWO-WAY 
ANOVA (TREATMENTS WITH BLOCKING). (DEFAULTED 
GROUP: TWOWAY 

SUBPROGRAMS: INPUT, SUMSQ, FMT 

DESCRIPTION : ANOVA DOES A ONE-WAY OR A TWO-WAY ANALYSIS OF 

VARIANCE DEPENDING ON THE VALUE OF THE OPTION PARAMETER. ANOVA 
WILL INTERACT WITH THE USER TO OBTAIN THE REQUIRED INFORMATION 
AND DATA. WHEN ENTERING DATA, SEPARATE THE DATA POINTS WITH AT 
LEAST ONE BLANK AND USE NO OTHER DELIMITERS . OUTPUT CONSISTS OF 
AN ANOVA TABLE, ESTIMATES OF THE OVERALL MEAN, AND ESTIMATES OF 
THE TREATMENT (AND, IF 0PTI0N=2, BLOCK) EFFECTS. 



AZRHOW 

SYNTAX: A3R VECTOR 

GROUP: SMOOTH 

SUBPROGRAM: 0NE3 

DESCRIPTION : AZR DOES REPEATED SMOOTHINGS OF RUNNING MEDIANS OF 3 
UNTIL THERE ARE NO CHANGES IN THE SMOOTHED SEQUENCE FROM ONE 
ITERATION TO THE NEXT. SEE ONEZHOW . 



AZRSRHOW 

SYNTAX: AZRSR VECTOR 

GROUP: SMOOTH 

SUBPROGRAMS: AZR, ONEZ, SPLIT 

DESCRIPTION : AZRSR DOES REPEATED SMOOTHINGS BY RUNNING MEDIANS 

OF 3 FOLLOWED BY SPLITTING MESAS AND REPEATING UNTIL CONVERGENCE. 
SEE AZR AND SPLIT. 
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BOXPLOTHOW 

SYNTAX: BOXPLOT VECTOR 

PARAMETER: 

LENGTH’ CONTROLS THE HORIZONTAL SIZE OF THE DISPLAY 
(.DEFAULT-^Q CHARACTERS). 

GROUP: DISPLAY 

SUBPROGRAM: FILL 

DESCRIPTION : BOXPLOT GENERATES A BOX PLOT DISPLAY FOR A VECTOR 

OF DATA. A RECTANGULAR BOX NITH ENDS CORRESPOND! NG TO LONER AND 
UPPER QUARTILES IS PRESENTED WITH THE MEDIAN MARKED WITH AN 
ASTERISK. NHISKERS ARE DRANN ON EACH SIDE OF THE BOX NITH 
CROSSES MARKING THE LOWEST AND HIGHEST DATA VALUES WIT til N AN 
INTERQUARTILE DISTANCE OF THE QUARTI LES . DATA VALUES OUTSIDE THE 
CROSSES {OUTLIERS) ARE MARKED WITH CIRCLES AND THOSE MORE THAN 
1.5 INTERQUARTILE DISTANCES GET HEAVY CIRCLES. 



CHISQU AREHOW 
SYNTAX: T CHISQU ARE X 

GROUP: GFIT 

DESCRIPTION: CHISQU ARE COMPARES A THEORETICAL DISCRETE 

PROBABI LITY MASS FUNCTION T HAVING N POSSIBLE VALUES WITH Ay 
EMPIRICAL FREQUENCY FUNCTION X OVER THE SAME N VALUES TO TEST HO^ 
WELL THE PROBABILITIES T FIT THE DATA X. THE Cil-SQUAtxE 
GOODNESS’OF’FIT STATISTIC IS PRINTED OUT WITH THE NUMBER OF 
DEGREES OF FREEDOM. 



CODERESHOW 

SYNTAX: CODERES MATRIX OR CODESRES MSDPOLISH MATRIX 

SUBPROGRAM: AND 

GROUP: TWOWAY 

DESCRIPTION: CODERES PRODUCES A DIAGNOSTIC ARRAY FUR A'NALYZIM', 

THE residuals FROM MEDIAN POLISHING OF A TWO’WAY TABLE. THE USt.'t 
SHOULD SET RESIDS-1 IN MEDPOLISii WHEN CODERES IS USED. Tue 
OUTPUT MATRIX CONSISTS OF THE SYMBOLS -°+ DEPENDING ON THE SIZE 
OF THE RESIDUALS. 

>: VALUE BELOW LOWER QUARTILE 

°: VALUE BETWEEN QUARTILES 

+: VALUE LARGER THAN UPPER QUARTILES 
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COMPAREHO^fi 

SYNTAX: COMPARE MATRIX 

PARAMETERS : 

(1) DEPTH- VERTICAL HEIGHT OF DISPLAY iDEFAULT=20 
LINES) . 

(2) MISS- NUMBER USED TO INDICATE MISSIVG VALUES 
(DEFAULT =~99999) 

(3) NGAP- NUMBER OF HORIZONTAL SPACES BETWEEN THE 
BOXPLOT DISPLAYS (DEFAULT=3) . 

GROUP: COMPARISONS 

SUBPROGRAM: FILL 

DESCRIPTION : COMPARE OPERATES ON AN N BY K MATRIX TO PRODUCE K 

VERTICAL BOX PLOTS PLACED NEXT TO EACH OTHER TO ALLOW VISUAL 
COMPARISON OF THE CENTERS, SPREADS, AND OUTLIERS OF THE BATCHES 
(COLUMNS) OF THE MATRIX. THE USER MUST FILL UP THE MATRIX SO 
THAT THERE ARE N OBSERVATIONS FOR EACH BATCH. THE PARAMETER, 

MISS, SHOULD BE USED TO FILL IN MISSING VALUES. THE USER SHOULD 
ASSURE THAT MISS IS DIFFERENT FROM ALL VALID ENTRIES IN THE DATA 
ARRAY . 



CONDENSEHOW 

SYNTAX: CONDENSE MATRIX OR R^CONDENSE MATRIX 

PARAMETERS : 

(1) NUM- CONTROLS WHAT STATISTICS ARE INCLUDED IN THE 
SUMMARY. /V£/AJeCl,2, 5,7] (DEFAULT = 2). 

NUM=1 GIVES MEDIANS OF EACH COLUMN OF THE MATRIX 
NUM=2 GIVES MEDIANS AND INTERQUARTILE RANGES 
NUM=S GIVES MIN,QUARTILES,MAX, AND SIZE 
NUM=1 GIVES MIN , EIGHTHS , MAX AND SIZE 

(2) MISS- NUMBER USED TO CODE MISSING VALUES 
( DEFAULT=~S9999) . 

GROUP: COMPARISONS 

SU BPROGRAMS: UTCOND, FMT 

DESCRIPTION : CONDENSE GENERATES SUMMARY STATISTICS FOR EACH 

GROUP OF DATA REPRESENTED BY THE COLUMNS OF THE ARGUMENT MATRIX. 
SINCE THE ARGUMENT IS A MATRIX, THE USER MUST FILL UP THE MISSING 
VALUES IN THE MATRIX WHENEVER THE GROUPS (COLUMNS) HAVE DIFFERENT 
NUMBERS OF OBSERV ATIONS . THE PARAMETER, MISS, SHOULD BE USED TO 
FILL THE MATRIX SINCE THE FUNCTION WILL RECOGNIZE THOSE VALUES /^5 
MISSING DATA AND WILL IGNORE THEM IN ALL CALCULATIONS . THE 

RESULT FROM CONDENSE CAN BE USED AS THE ARGUMENT OF OTHER 
FUNCTIONS SUCH AS SCAT OR LINE IF AND ONLY IF NUM=2. 
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CONTINGENCYHOy 
SYNTAX: CONTINGENCY MATRIX 

GROUP: T'NOWAY 

DESCRIPTION : CONTINGENCY TAKES A TNO-WAY TABLE AND PERFORMS A 

TEST OF INDEPENDENCE OF ROWS AND COLUMNS. THE CHI~SU'JARn 
STATISTIC AND ITS DEGREES OF FREEDOM ARE OUTPUT. 



KSHOW 

SYNTAX: F KS G 

GROUP: GFIT 

DESCRIPTION : KS DOES A KO LMOGOROV ~SMI RNOV GOODNESS OF FIT TEST 

OF THE EQUALITY OF TWO DISCRETE EMPIRICAL PROBABILITY 
DISTRIBUTIONS . THE ARGUMENTS F AND G ARE EACH VECTORS OF 
OBSERVATIONS. THEY NEED NOT BE SORTED, NOR MUST THEY BE THE SAME 
SIZE. OUTPUT OF KS ARE THE VALUES: 

MAX F(X)’GiX) AND MIN FU)-GiX) 



LINEHOW 

SYNTAX: LINE ARRAY OR Z^LINE ARRAY 

PARAMETER: 

RSELECT’ SELECTS THE OUTPUT MATRIX OF RESIDUALS. 
RSELECT-1 GIVES ABSCISSA VS. RESIDUALS {DEFAULTED . 
RSELECT*! GIVES FITTED VALUES VS. RESIDUALS 
GROUP: RELATIONS 

DESCRIPTION : LINE FITS A STRAIGHT LINE TO A SET OF iX,Y) POINTS 

BY DIVIDING THE POINTS INTO 3 REGIONS AND USING THE MEDIANS OF 
THE X AND Y VALUES IN THE OUTER REGIONS TO DETERMINE THE SLOPE. 
THE INTERCEPT IS THE MEDIAN OF THE DIFFERENCES Y^SLOPE^X. Llut 
ALSO COMPUTES THE RESIDUALS AND GIVES AN N BY 2 MATRIX CONTAINING 
EITHER THE X-VALUES VS. RESIDUALS OR THE FITTED VALUES VS. 
RESIDUALS DEPENDING ON THE PARAMETER, RSELECT. IF THE USER DOES 
NOT WANT RESIDUALS TYPED OUT, HE SHOULD USE THE SYNTAX: 

Z^LINE ARRAY 

THE ARGUMENT ARRAY CAN BE EITHER A VECTOR OR A 2~C0LUMN MATRIX. 

IF THE ARGUMENT IS A VECTOR OF SIZE N, THE X^VARIABLE IS 
CONSTRUCTED TO BE THE FIRST N POSITIVE INTEGERS. IF THE ARGUMENT 
IS A 2-COLUMN MATRIX, THE FIRST COLUMN IS TAKEN TO BE THE SET OF 
X-VALUES AND THE SECOND COLUMN THE SET OF Y-VALUES. LINE CAN BE 
USED IN CONJUNCTION WITH THE FUNCTION SCAT TO PRODUCE THE SLOPE, 
THE INTERCEPT, AND A PLOT OF THE RESIDUALS. TO DO THIS, ENTER: 

SCAT LINE ARRAY 
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LSLINEHO^ 

SYNTAX: LSLINE ARRAY OR Z^LSLINE ARRAY 

GROUP: RELATIONS 

DESCRIPTION : LSLINE DETERMINES THE LEAST SQUARES SOLUTION, 

Y-YI*SLOPE^X . THE ARGUMENT ARRAY CAN BE A VECTOR OR A 2^C0LUMN 
MATRIX. IF THE ARGUMENT IS A VECTOR, THE X’VALUES ARE TAKEN TO 
BE THE INTEGERS 1 TO N WHERE N^p ARRAY. IF THE ARGUMENT IS A 
2~C0LUMN MATRIX, COLUMN 1 CONSISTS OF THE X^VALUES AND COLUMN 2 
THE Y^VALUES. RESIDUALS ARE DETERMINED AND STORED IN A MATRIX 
WHOSE FIRST COLUMN IS THE SET OF X ’VALUES AND WHOSE SECOND COLUMN 
IS THE SET OF RESIDUALS. IF THE USER WANTS TO SUPPRESS PRINTOUT 
OF THE RESIDUALS, HE MUST ASSIGN THE RESULTS OF LSLINE TO A 
VARIABLE; I.E., TYPE Z^LSLINE ARRAY. 



MEDPOLISHHOW 

SYNTAX: MEDPOLISH MATRIX OR Z^MEDPOLISH MATRIX 

PARAMETERS: 

( 1 ) RESIDS’ CONTROLS OUTPUT OF RESIDUALS. IF 1, 

RESIDUALS ARE PRODUCED IN A TWO-WAY TABLE; IF 2, 

COMPARISON VALUES, iRE° .^CE) iTV AND RESIDUALS ARE 
PRODUCED IN A TWO COLUMN MATRIX (DEFAULTED . 

( 2 ) EPSILON’ PROPORTION BY WHICH THE SUM OF THE 
ABSOLUTE VALUES OF THE RESIDUALS MUST BE REDUCED AT 
EACH ITERATION TO CONTINUE POLISHING { DE FAU LT = 0 . 0 1) . 

( 3 ) NORMEFFECTS’ DETERMINES WHETHER OR NOT NORMALIZED 
EFFECTS ARE OUTPUT. THE DEFAULT OF 0 SUPPRESSES OUTPUT 
OF NORMALIZED VALUES. N0RMEFFECTS=1 CAUSES NORMALIZED 
VALUES TO BE PRINTED. 

GROUP: TWOWAY 

DESCRIPTION : MEDPOLISH ITERATIVELY SWEEPS OUT MEDIANS FROM THE 

ROWS AND COLUMNS OF A TWO-WAY TABLE TO YIELD THE MODEL: 

OBS = MEDIAN + ROW EFFECT + COLUMN EFFECT + RESIDUAL 
IT YIELDS RESIDUALS FOR TESTING THE ADEQUACY OF THE MODEL (SEE 
SHOWRES AND CODERES) . THE FUNCTION CONTINUES POLISHING UNTIL THE 
STOPPING RULE CONTROLLED BY EPSILON IS ACTIVATED. THE OVERALL 
MEDIAN, THE ROW- AND COLUMN-EFFECTS, AND THE SUM OF ABSOLUTE 
VALUES OF THE RESIDUALS ARE GIVEN. THE USER CAN SUPPRESS 
RESIDUAL PRINTOUT BY ASSIGNING THE FUNCTION TO A VARIABLE, I.E., 
BY WRITING: 

Z*- MEDPOLISH MATRIX 
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CORRELATIONHO^</ 

SYl^TAXt R^CORRELATION Iv' 

GROUP i STATS 
SUBPROGRAMt FMT 

DESCRIPTICti \ CORRELATION DETERMINES THE SIMPLE PEARSON PRODUCT 
MOMENT CORRELATIONS BETWEEN EACH PAIR OF VARIABLES REPRFSt NTED BY 
THE C COLUMNS OF W, THE OUTPUT IS A C BY C CORRELATION MATRIX. 



MSTATSHOW 
SYNTAX t MSTATS W 
GRCUPi STATS 
SUbPhOGRAMi FMT 

DESCRIPTION \ MSTATS DETERMINES CORRELATIONS , MEANS, STANDARD 
DEV lATIONS , LOWER QUARTILES, MEDIANS, AND UPPER QUARTILES FOR THE 
COLUMNS OF THE ARGUMENT MATRIX W. THE CORRELATIONS ARE DISPLAYED 
IN AN N BY N MATRIX, WHERE N IS THE NUMBER OF COLUMNS OF W. 



PARTI ALHOW 
SYNTAXi PARTIAL W 
GROUP t STATS 

SUBPROGRAMS t CORRELATION , FMT 

DESCRIPTION \ PARTIAL DETERMINES PARTIAL CORRELATIONS OF ANY 
SPECIFIED ORDER FOR THE COLUMNS OF THE MATRIX W. THE USER WILL BE 
ASKED FOR THE COLUMNS TO BE PARTIALLED OUT. THE USER'S RESPONSES 
MUST BE INTEGERS BETWEEN 1 AND C, WHERE C IS THE NUMBER OF 
COLUMNS OF W. 
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NUMSUMHOW 

SYNTAX: NUMSUM VECTOR OR NUMSUM ARRAY 

PARAMETER: 

MISS’- NUMBER USED TO INDICATE MISSING VALUE IN THE DATA 
ARRAY ( = "99999 ) . 

GROUP: COMPARISONS 

SUBPROGRAM: FMT 

DESCRIPTION : NUMSUM OPERATES ON EITHER A VECTOR OR AN ARRAY BUT, 

IN ALL CASES, THE DATA ARE TREATED AS A SINGLE BATCH. NUMSUM 
PRODUCES A NUMERICAL SUMMARY WHICH GIVES THE SAMPLE SIZE {AFTER 
DELETION OF MISSING VALUES), THE LIGHTHS, THE EXTREMES, THE 
SPREADS, AND THE MIDPOINTS IN TABULAR FORM. THIS SUMMARY IS 
USEFUL FOR TESTING THE SYMMETRY OF A DATA SET, AND TO EVALUATE 
THE EFFECTIVENESS OF DATA TRANSFORMATIONS IN PRODUCING SYMMETRY. 



ONEHHOW 

SYNTAX: ONEH VECTOR 

GROUP: SMOOTH 

DESCRIPTION : ONEH DOES A SINGLE HANNING OF A SEQUENCE OF DATA. 

IT IS USED IN CONJUNCTION WITH THE OTHER SMOOTHING FUNCTIONS AS A 
FINAL TOUCH-UP TO A SEQUENCE OF SMOOTHED DATA. THE I-TH RESPONSE 
IS REPLACED BY 0.25xY[J-l] + 0.5xy[J] + 0.25xy[J+l]. SEE AZR AND 
A2RSR. 



0NE3H0W 

SYNTAX: 0NE3 VECTOR 

GROUP: SMOOTH 

DESCRIPTION : 0NE3 SMOOTHS A VECTOR OF DATA USING ONE PASS OF 

RUNNING MEDIANS OF 3 WITH TUKEY'S END-POINT RULE. WHEN SMOOTHING 
(X,Y) PAIRS, THE USER SHOULD USE AS THE ARGUMENT VECTOR THE Y 
VALUES ORDERED ACCORDING TO THE MAGNITUDE OF THE X VALUES: 

VECTOR^YliX^ 

IT IS ASSUMED THAT THE X VALUES ARE EQUISPACED. REPEATED 

SMOOTHING BY MEDIANS OF 3 CAN BE ACCOMPLISHED WITH A3R. OUTPUT 
CONSISTS OF THE SMOOTHED SEQUENCE. A PLOT OF THE SMOOTHED 
SEQUENCE IS GIVEN BY TYPING: 

SCAT ONEZ VECTOR 
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REGREFSHOW 

SYNTAXx Z^Y REGRESS X 
P ARAMETERx 

ElNTERCEPTr. DETERMINES WHETHER OR HOT AN INTERCEPT Tt.RM 
IS TO BE INCLUDED, hlNTERCEPT^l GIVES AH INTERCEPT 
TERM, AND hINTERCEPT-0 GIVES NO INTERCEPT, {DEFAULT IS 
1 .) 

GROUPi RELATIONS 
SUBPROGRAMSi FMT AND SCAT 

DESCRIPTION t REGRESS DOES A MULTIPLE REGRESSION ANALYSIS 
RELATING THE DEPENDENT VARIABLE Y TO A SET OF CARRIERS X. THE 
LEFT ARGUMENT Y IS A VECTOR OF SIZE N, THE RIGHT ARGUMENT X IS 
AN N BY K MATRIX CONSISTING OF N OBSERV ATIONS ON EACH OF K 
VARIABLES OR A VECTOR OF SIZE N IF /C = l. OUTPUT CONSISTS OF /1,-v 
ANOVA TABLE, RmSQUARE, STD, ERROR, REGRESSION COEFFICIENTS {Tdt 
FIRST COEFFICIENT IS THE CONSTANT TERM IF M NT E RC E PT= 1 , ) , 

TnSTATISTICS , V ARI ANC E r,C OV ARI ANCE MATRIX, DURBIN ATSON 

STATISTIC, AND A VECTOR OF PREDICTED Y VALUES AND RFSIDUALS, 
THERE IS AN OPTION THAT ALLOWS THE USER TO INPUT A VECTOR OF X 
VALUES AND USE THE REGRESSION EQUATION TO FORECAST Y VALUES, THE 
USER CAN ALSO OBTAIN A SCATTER PLOT OF THE RESIDUALS, WltFi'i 
EXECUTION TERMINATES , THE PREDICTED Y VALUES AND THE RESIDUALS 
RESIDE IN THE N BY 2 MATRIX Z, 



STATISTICSHOW 

SYNTAXx STATISTICS VECTOR 
GROUP X STATS 
SUBPROGRAMx FMT 

DESCRIPTION X STATISTICS DETERMINES THE MEAN, VARIANCE, STANDARD 
DEV I ATION , COEFFICIENT OF V ARIATION , LOWER AND UPPER QUARTILES, 
MEDIAN, TRIMEAN, MIDMEAN, MIDRANGE, RANGE, MEAN ABSOLUT R 
DEVIATION, INTERQUARTILE RANGE, SKEWNESS, AND KURTOSIS EOR THt. 
DATA IN VECTOR, 
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SCATHCi>i 

SYNTAXi SCAT ARRAY 
PARAMETERS X 



(1) WID-^ CONTROLS THE HORIZONTAL SIZE 
{DEFAULT-ZO CHARACTERS) , 


OF 


THE 


DISPLAY 


(2) DEPr CONTROLS THE VERTICAL SIZE 
{DEFAULT=\b LINES), 


OF 


THE 


DISPLAY 


(3) NDIVX, NDIVYv NUMBER OF UNITS ON 
RESPECTIVELY ( D5F-4i'Lr = 4 , 4 ) . 

GROUP: RELATIONS, SMOOTH, COMPARISONS 


Xr- 


AND 


Y ’■AXES, 



DESCRIPTIONx SCAT PRODUCES A SCATTER PLOT OF THE DATA CCNT AIRED 
IN ARRAY, THE ARGUMENT ARRAY CAN BE A VECTOR OR A MATRIX WITH -45 
MANY AS 9 COLUMNS, IF A VECTOR OF SIZE N, THOSE VALUES ARE 
PLOTTED VS, THE INTEGERS 1 TO N\ IF A MATRIX, THE SECOND, THIRD, 
ETC, COLUMNS ARE PLOTTED K5. COLUMN 1 ON THE SAME AXES, 
DIVISIONS ON THE AXES OF THE PLOT ARE NOT EXPLICITLY PRINTED, 
EXCEPT AT THE EXTREMES OF THE PLOT, THE USER CAN CONTROL THE 
RESOLUTION OF THE PLOT BY MODIFYING THE PARAMETERS WID, DEP, 
NDIVX AND NDIVY, PRINTING TIME INCREASES DRAMATICALLY WITH 
RESOLUTION -, THEREFORE, OTHER PLOT PROGRAMS SHOULD BE USED IF HIGH 
RESOLUTION IS DESIRED, FOR A SINGLE GROUP OF DATA, THE NUMBER N 
i2<N^S) WILL BE PRINTED IF N POINTS LIE CLOSE TOGETHER ON THE 
DISPLAY, FOR MULTIPLE PLOTS ON THE SAME DISPLAY, THE LETTER A 
REPRESENTS GROUP 1, B REPRESENTS GROUP 2, ETC, FOR 2 POINTS 
CLOSE TOGETHER THE LETTER WILL BE PRINTED WITH AN UNDERSCORE , 
THE DISPLAY CANNOT HANDLE 3 OR MORE POINTS CLOSE TOGETHER IF 
THERE ARE MULTIPLE PLOTS, 
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SHOWRBSHO'^l 

SYRTAK: SHOWRBS MATRIX OR SHOWRES MEDPOLISH MATRIX 

SUBPROGRAM'. AND 
GROUP: TWONAY 

DESCRIPTION : SHOW RES AIDS IN THE ANALYSIS OF RESIDUALS FROM 

MEDIAN POLISHING A TWO-WAY TABLE. THE USER SHOULD SET RESIDS=1 
IN MEDPOLISH IF SHOW RES IS TO BE USED. THE OUTPUT OF SHOW RES IS 
A MATRIX CONSISTING OF THE SYMBOLS OOoxX DEPENDING ON THE SIZES 
OF THE RESIDUALS. LET R REPRESENT THE RESIDUAL, M THE MIDSPREAD 
OF ALL RESIDUALS, L THE LOWER QUARTILE, AND U THE UPPER QUARTILE. 
THE CODES ARE PRINTED AS FOLLOWS: 

0: R<L~1.SM 
O; L~1.5MiR<L^M 
L’M<R<U+M 
x: U+M<R<U+1.SM 
X: R>U+1.5M 



SPLITHOW 

SYNTAX: SPLIT A3R VECTOR 

GROUP: SMOOTH 

DESCRIPTION : SPLIT DOES ONE PASS AT DIVIDING MESAS {PAIRS OF 

ADJACENT POINTS WITH A COMMON VALUE WHICH IS A LOCAL MAX OR MIN) 
USING TUKEY'S END-POINT RULE. IT IS USED IN CONJUNCTION WITH 
A3R. 
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STSMLBAFHOW 

SYtITAX: STEMLEAF VECTOR 

GROUP: DISPLAY 

PARAMETERS : 

( 1 ) WIDTH’ CONTROLS THE WIDTH (CHARACTERS PER LIVE) OF 
THE DISPLAY ( DEFAU LT=7 0) . 

( 2 ) SCALE’ VARIES THE DEPTH (STEM INTERVAL) OF THE 
DISPLAY IN UNITS OF 1 , 2 , OR 0.5 TIMES A POWER OF 10 . 
SELECT AN INTEGER FROM 1 TO Z (DEFAULTED . 

SU tiPROGRAM: LIT 

DESCRIPTION : STEMLEAF GENERATES A STEM AND LEAF DISPLAY OF A 

VECTOR OF OBSERV ATIONS . THE FUNCTION AUTOMATICALLY SCALES THE 
DATA USING A SCALING ROUTINE BASED ON THE RANGE AND SIZE OF THE 
DATA VECTOR. BY VARYING THE SCALE PARAMETER THE USER CAN CHANGE 
THE SIZE OF THE STEM, BUT IT WILL NOT NECESSARILY BE SCALED BY 
THE AMOUNT SPECIFIED. TRY VARIOUS VALUES LIKE 1 , 2 , AND 3 TO SEE 
WHICH CHOICE YIELDS THE BEST DISPLAY. IF A LEAF CONTAINS TOO 
MANY CHARACTERS TO BE PRINTED ON A LINE, THE LINE WILL BE 
TRUNCATED AND THE NUMBER OF TRUNCATED CHARACTERS WILL BE WRITTEN 
AT THE END OF THE LINE. IF THE STEM INTERVAL IS TWO TIMES A 
POWER OF TEN, THE LEAVES MAY CONTAIN THE CHARACTERS A,B,C , . . . ,-J 
IN ADDITION TO THE DIGITS 0 , 1 ,..., 9 . IN THIS CASE, THE 
CHARACTERS A TO J REPRESENT THE NUMBERS 10 TO 19 . 



TWICEHOW 

SYNTAX: TWICE VECTOR 

GROUP: SMOOTH 

SUBPROGRAMS: ONE3, SPLIT, A3RSR 

DESCRIPTION : TWICE SMOOTHS A SEQUENCE OF DATA, THEN SMOOTHS THE 

RESIDUALS AND ADDS THE SMOOTHED RESIDUALS BACK TO THE SMOOTHED 
DATA TO OBTAIN THE FINAL SMOOTH. THE SMOOTHING IS DONE USING 
A3RSR. 
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9>VEFAULT values of TilE PARAMETERS 





OFF 


1 5 


DEPTH 


20 


EPSILON 


0.01 


LENGTH 


50 


MISS 


"99999 




tvDIVX 


4 


NDIVY 


4 


UGAP 


3 


iWRMEFFECT 


0 


HUM 


2 


OPTION 


1 


KESIDS 


1 


SCALE 


1 


RSELECT 


1 


NID 


30 


NIuTH 


70 


MNTERCEPT 


1 


CH 
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APPENDIX B 
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^AiwL'd:\^ 

V L^A AMD B;C;D 

Cl] ■>( ( (2<pp/i)v3<pp£^),0;^ppB)/ 17 3 

C 2 ] B -*- , b 

[3] ->( ( (3 = pp£()Ai?:ippB) ,2 = pp/^ )/ 17 7 

[4] A<-^A 

[5] ->(A/((p/^)^l,0),l^i)-*-lp“24)p6)/16 

[6] ( (Z^xp^)L/;rp/J) ,l)p/l 

[7] -y(lxppt) / 9 

[8] fl-^( ( ( p6 ) r ( 1 =p^ ) X 1 pp/1 ) , 1 ) plJ 

[9] ->( (A/Z?;^l,lpp^),lxZ)^lp~24>P«)/ 16 11 

[10] (.(3 = ppB)pl)^ilppA)ylpi>pB)pB 
Cll] ->(3 = pp6)/l4 

[12] L-^(((r-<-lpij)p/4)p0),(lp<t>pii)pl)\fl 

[13] ->0xppi[ ; 

[14] ( 1 , ( ( T"*-! p 4>P/4 ) p 0 ) , ( 1 + 1 p<t5p6 ) p 1 ) \fi 

[15] -+0 xp pL[ ; ; 1 + iC ]-^^ 

[16] ^0 = pLl^' ARGU'uBMTS OF AMD AkB MOT COMFORMABLf. ' 
L17] 'AN ARGUMBMT OF AMD IS OF IMRRORBR RANK.' 

V 



[ 1 ] 

[ 2 ] 



MAMOVAi:\']M 



V AMOV A \Y \K\BTSS \‘MSR\R\MAK\ BBSS \yiSBL'J'lSB‘,MT S3 \S OF \TDF\T SB -,1 BA 



Ti BLOCK 



INPUT 



SUMSU 

V 



V/!|3i^[a]V 

V Z^A3R X\Yl;Y2;Xl;X2iM-,C 

[1] Z->-OMS3 X 

[2] -vix ( ppp^^Z) x0^ + / I 

V 



MA3RSRLU1M 

V Z^A3RSR X 

[1] X-^A3R X 

[2] Z^A3R SPLIT X 

[3] ->-2x(pppA'-<-Z)xO== + /|Z-’A' 

V 
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vzjoa:ploj’[0]v 

V Z^bOXPLOT X-,L-,U‘,b 

[ 1 ] X-i ill!:~20+U ^L)^{ iU-^{ / ,X)^Le!VGTii^L-^l / ,X ) + ^X^LbNGTH^l 

[2] 0. IxLO. 5 + lOxZ/; ( orLfiA/GTrt-e )p ' ' ; 0 . 1 x [ 0 . 5 + 1 Qxi/ 

[3] -*-9xxO=y’L 

[4] Z^LENGTb FILL X 

[53 P-' "I n r ’[1+Z] 

[6] B-^b,' *1*1* >x ooooooooo ' [ i + Z ] 

[7] _ILII_ 23456789 2 3 4 5 6 7 8 9 ' [ 1 +Z ] 

[8] -►0xppZ^(3,Li’/l/G?’yy)pS 

[9] 'RAi\IGE OF DISTRIBUTION IS ZERO' 

V 



VCiilSQUAREiaiV 

V T CRISWARE X-,DF\C 

[1] -*-3x1 ( (p^)=py) 

[23 LENGTH ERROR: X AND T MUST BE THE SAME SIZE.' 

[33 C^^niX-‘T)*2)iT^T^*/X 

[43 'CHI--SQUARE VALUE = ';C;' DF - ';(pZ)’l 

V 



NCODERESlW^N 

V Z->^CODERES X ;Y ;N \Q;LN iUU 

[13 /i?-py^y[^y^,i3 

[23 3-^0.5xy[L0.2 5x3+/Vxi33+y[r 0 . 25 x 1 +^ 7 x 13 ] 

[33 [ 1 3 [ 3 3 [ 1 3 

[43 UN^QLz:i+D 

[53 y^{X>LN)-^{X>QLll) + (X>QL3'])+X>UW 
[63 Z^’’»o++' [i+y3 

[73 Z^((pZ[;l3),2xpZ[l;3)p(,Z) ANDApZ)p' ’ 

V 



94 



'^COMPARt'LQly 

7 Z^COMPARE X W iliy iJ iP 

Cn X^iilE 20-^i'-Q-^L) x( /Y)-DFPTH^L-^l/Y*-U'}ISS* ,X) / ,X)-rX^ht.P':n 

[ 2 ] W-e- ( T If ■ 2 0-ri' -L) X ( 6' nDEPTH^L ) -^-MISS^DEPTH-^l 

[31 Z-^(DfP'i'//,6xC’^pA'[l; Dp' ’ 

[4] J-^0 

[ 5 1 J '*-J + 1 

[6] P^DEPTH FILL PLiP^(P^M) /P^XL -,J ']} 

[7] ZC;"2-r6x^>’ ri_T| '[l^Pl 

[8l Z[ ; ~l + 6 x«7 1-*- ' *T = i= lx ooooooooo ' [ l -rP ] 

[9] Z[;6xJl-e« I"i_x| 23456789 2 3 4 5 6 7 8 9 ’ C 1 -rP 1 

[10] 

[111 Z-^«Z 

7 



7r(7/VPf/VP£'[G17 

7 R^C ON DENSE X \1 \N \N R ',U iQ \Z 
[ll P^lO 
[2l UTCCND X 
[31 -*0-<\{NUM*2) 

[4l R^{NR,2)qZ 
7 



VCONTINGENCYiUlSl 
7 CONTINGENCY F ;R -,C iT ;E ‘,C HI2 ;DF 
[ll R->- + ^F 

[2] T^+IC^+IF 

[3] E^Co,xRiT 

[4l CHI2^-^/ ,((FnE)*2)iE 

[51 PF^x/(pP)-l 

[6l TPI-PQtMPf l//1LPP= ';CHI2;' DF= ’;PP 
7 
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SICCNDENSEIU']'^ 

V R-<rCOt} DENSE X ;I ;N ;N H iU ;Q iZ 
Cn R^\0 
C 2 ] UTCOND X 

[3] ->~Q^\{NUMx2) 

[4] R^iNR,2)pZ 
SI 



SICONTINGENCYLaiV 
SI CONTINGENCY F ‘,R ;C iT ‘,E ‘,C HI2 DF 
Cl] R^+JiF 

[ 2 ] T^+/C^+/F 

[3] E‘*-Co,xRiT 

[4] CHI2^->-/ ,((F”E)*2) iE 

[ 5 1 DF-*-^/ ( pF) r 1 

[6] 'CHImSQUARE VALUE= 'iCHI2;' DF= *;£»F 
7 



^CORRELATIONiniV 

7 R^CORRELATION W ;Z iC ;S iCgiMEANS iV AR 
Cl] C^{isiZ)'^.^Z^Nr,(pW)p{MEA'NS^ + ^W)il^pi~2fl,l,QW)pW->-pCH^' ' 

C2] R^'BFS,2' FMTiO,\pS) Al'UspS) ,CiSo,xS->-(VARS^->-/,Z*2)* 

0.5 
7 
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'JFMTIQIV 

7 OL^l FMT RiS\W \L\G\X\T \K\J \0\L\B\V \CH\H 

[1] 2\qR)qR 

[2] = 1 0 + ' , ' )pA-^’0123456'’89. ' 

[3] ->£’xi (/V^O=.V) vl/-Hi>pS-H,£ 

[4] t0:-^~i7v( xp^4x2 = p/^^pjf^» ' ) Av/( ' ,(9^»D' )eS 

[5] ->(Zy0 + (F-<-0=p5-<-J^ + 5)+xS=;W[2]-<-l),_^-( \'x-D -*-0 + . -K) A ' eK*-K , ( •/ -^5 x 

+ 5 

[6] -yE^^-QS^'TEXT DELIMITER^ 

[7] -►Z/ 3 » 3 X X ( p(j -*-/(=/(;■<-( ^ £ Hr E) / K ) \.W-*-p X-*-( p K-*-{ K \ 0 ) i K ) ^ i ” { i) K ) \ 0 ) ^ K 

[8] L:^{D^~HG->-K€A) /L3’2^ipK)*W^lW^' XA' €K^{~Ke' H)/K 

[9] ->L3x t (p?: + /G )«xM[ 2]-^-10i I l>Ax ( 1>(7 X 0) f A' 

CIO] ->A3-<t)0, -(L-^-’AA/’eA ) / xA/^-lOl 1 1-Ax ( | l-G-<-£ix H ' ) i B^(. ( i>G ) \ 0) i K 

C 11] l + pA-^-C ( Ifpi) L (MC 1] >//) , V) 

Cl2] L3:->(A£ixxAA~»A'eA),e-pA-^-l7,D-^0pP'^( (/7-tf , 0 ) x i , V ) pA 

[13] -^A4-x~ltA,S'^l+pP^(0 1 xpp^pc ; xA/[2]-^C)L>iC2]rCx[/AO])4-A 

[14] P«^Pt10*L>^[ 10® I P+0 = P 

[15] ->A3xx 0=./^ + /P^( ’ S' eA)x>t0=P-^( [0 . 5+,Vx,p) i/i/^10*P-^10i I 1-A xG + 0 

[16] A4:->(pl4-pA)/P-ppA>^(l 0 ^pG*-JpT\' ' ) pJ *-J ,0*-^ /T*-0>P^B / P 

[17] •^-(xp-K(orAxj.^»Z’eA)[. x-y-eC y + 0.^-1 + L 1 0® 1 f | P ) >0 + -P +0 +~2 i A ) / A/ 

I 

[18] -^E + ^-pS^' FIELD WIDTH' 

[19] ->A4 + l + x ( (J[ 2]-^Av. <0)+G^l + 10f . < | L*-i B / , L) +T^10= \ P ) >W ’D +0 + 

3 

[20] T-^'^J*-P[T / X lf-7 ]-^A'^plpA'^-'P’ , ' + 0~ * [•/ p 2 - xp ] , A [ 1 + 5? ( Op 1 0 ) T i A ] 

[21] P:->(Jv2>Ax~*7’»eA)/J,A/^pA^A[ll,l+fi?(AplO)Tr/Zxl| |P],A 

[2 2] P-^, ( -P) i ( A[ . xSjaC ; 2+A]xl+A ) o . <P-^xP»l 

[23] A^/l/pA,A[P/xpA-^,A]'>-' ' 

[24] J vO = + /G-^OrA-P)/A + pP-^pP-^G, A[1 + <S?(Ap10)tL | P] 

[2 5] P^Pp ( ,0^G4>P)\( ,P^po.<(-G)4>xA+G'^l^pG)/,P 

[2 6] ->PP-xJ'vA^~' A' eA,P[J'/ xP^ltA^pP^P, A; ]^' * ' 

[27] P^Ap ( ,<t>P)\( ,P-^~A + ~0)/,P 

[2 8] •>(~P)/A-/i/^l,P^OpP>^A\(P, A-^Vxl-2xA)+P 

[2 9] AP :rA-^( pA^( "l+P-^O, (A; [ 2 ] [ pP ) p P-^ ( ’ , • =C H ) / \ pC H ) pC H )<i>C ,i 

[30] P-^-, (AZ[2],A)t 0 ~1 ^(A;[2],e)p(,<tP‘>.>xA-^-r/P“^l + P-~l<t>P)\A 

[31] ->(A0-1/ax3) , pGPhkGP, ( (l = l+A'mA/xi,^)pP, ,P 

[32] E:K^'H0 VALID E, 7, OR F PHRASE' 

[33] 'FMT PROBLEM ',A;(l,pS)pS 

7 
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7FIZ, LCDIV 

7 FILL XiN iDiXUX2;X3iXHiY il iJ iS 

Cn U^pX^XL^X^.X'l 

[2l Z■^0.5x;'CL0.2 5x3+^x^3]+;^[^ 0 . 2 5 x 1 +.V x i 3 1 

[31 5^(Z[ll-0. 5 + 1. 5xZ?) . (Z[n^0.5+Z?-^Z[ 31 nZ[ n ) .ZCll 

[4l 5'+L0.5+5.ZC 21 . ZC3l.(ZC3l + 0.5+D),Z[3l■^0.5 + 

1 . 5 x£) 

[51 Xl-^(X2<SLll) /X2->-(X^SL2'\) /X-^10,5+X 

[6i A’2>-(;i:2>5[ 11 ) /;if2 

[71 a: 4-^U3>5[71 )/;f3-^U>5[6i )/;i; 

[81 ;f3-<-(;f3<5[7i )/;f3 

[9l y-+( (;f>5[ 21 ) xji:<5[6i ) / a: 

[101 Z-^VpO 

[111 Z[I>2 0 + 9L+,<("l,A'l) O.=l^\(f^0r5l;7) [11 

[I2I Z[lKlO-^9L+/("l,y2) o,=I'+5[ll + il+5[2lr,5[ll 

[131 Z[I1^10 + 9L + ,4Cl,y3) o .=I-^"l+5[6l + il+5[7l-»5[6l 

[141 Z[Il+-20 + 9L-^/y4o.=I-^.■■l^5[7l + ^l..V„5[7l 

[151 +18x\0=py 

[161 z[j-^"i+y]^z[i^i+y]^9 

[17] Z[I-mO[JtI + 11^8 

[18] Z[5[3l + \0r5[5l«.5[3l + ll*+6 

[191 Z[5[3l,5[5ll^0 

[20I Z[5[4ll+-1 

[211 Z[5[3].5[5]]^Z[5[3l ,5[5ll+ 2 4 

V 



7I/7Pt'r[Ql7 

V INPUT iJ iM;Z 

[1] 'ENTER NUMBER OF TREATMENTS, ' 

[ 2 ] J^0^K*-Q 

[31 'ENTER tKQ.ZQ.E. OF NUMBER OF DBS, FOR TREATMENTS 1 TO ' iK 
[4] y-( (M/5y^r7/V-D) ,^)p0 

[51 ->-7^\{K=pN) 

[61 ->-^yppQ->-' RE-ENTER COUNT VECTOR, ONE ELEMENT PER TREATMENT ,' 

[71 -►Ox I (/^<J^J + 1) 

[81 'ENTER 'iNLJli' OBSERVATIONS FOR TREATMENT ' -,J 
[91 +-7xpppy[ J 1^(MAX, l)pZ,(M^M AX r,pZ^Q)pO 

V 
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V 



E KE Li i E, i i\l ; X i l\i t' i N G i M AX ; f I u :v i 
[ 1 1 "-*-0 X •i’->-pZ-^Z [ ^i2+-p , J 1 

[ 2 l ->-j 4 1 X ; ( ;V < .-*-,.-rl ) 

C 3 1 yVf P , ( ^/F<x*-Z [ 1 ) f pP 

i.U] F , {f / 1 <X) i qF 

l 51 //CrX-(i-/Cr<A)fpcJ 

[ 6 1 FIGl\/->-y-,vF-<-,PjFo , - AG , ( ->-/:;<A' ) i pG 

C 7 1 0*-D , JiAX X 1 f ( ( i]/F=t',AX*-f / iiF-^ [ ) /PIG v ) 

L 8 1 -►2 , / V F-*- \ 0 

[ 9 ] Al i ' El AX F ( A ) - ;; ( X ) = ' ; r / A 
[JO] 'hlG PU)-G(X) =';L/G 
V 



7 LIAP[Q ]7 

7 INLINE W ;X ;X iN ;M;P;I iJ iK;Xl;YliX2;Y2iF;SLiYl 
Cll ^3xi2=ppV 
[2] I7^<s?(2,pl7)p(xpy),f7 

[31 ; n l)pA-^5L^0 

[41 X1^0.5xI[I^L0. 5x/>j>ll-rI[J^rO. Sxl-r/^-^-LO. 6^/0 31 

[5] Z2X-0. 5xI[A-rl-Il+I[A + lnJl 

[61 Yl-^ 0 . 5 xZ[Il + (Z-^Z[^Z-^M + y-^V[P; 2 lrFl)[Jl 

[71 y2^0.5xZ[lXZ-^-Z[U^( -A* )+Yl )[J1 

[81 F^(SL^SL^(Y2^Y1) iX2”Xl)^X 

[91 K->-K-rl 
[101 -► 6 xiA <2 

[111 yi^0.5xY[ L0.5xA + ll + (Yx-y[^Y^V[P;2lnPl) [fO.Sx/vTll 

[121 'SLOPEx ';5L;' Y ’^INTERCEPT i ';YI 
[131 Z-^y[ ; 2l-P-^-yi-rPLxfV[ ;il 

[141 -►lex iA5£'LiTP = l 
[151 ^0. ,Z^4?(2,A)pP,Z 

[161 Z^i)(2.A)pA[;ll.Z 

7 



7Lir[ai7 

7 R^LIT A 

[11 P^,' 0123456 7 89* [l-r( ( 1 + [ 1 0®/? ) p 1 0 ) T ^ 1 

7 
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VLSLJ/V£’CG] V 

V Z^LSLINE 
Cl] ->/^lxi2 = ppI 

[2] X^>S)i2^N)pX^{ iB^pX) ,X 

[3] ^l;6-;r[;2]@(6)(2»/V)p( ,XC; 1]) , (ly^lipX)pl) 

[4] 'SLOPE-. ’;B[1];' X^INT ERCEPT i ’;S[2] 

[5] 2^6!(2,:V)p(;r[;l]),^C;2]-.(B[2]+.?[ ;l]xe[l]) 

V 



IMEDPOLIStilU'I'y 

V Z-^MEDPOLISH X k ;C N ; RI -,C I \ R A \C A -,S \ RG RE -,C E -,TV 
Cl] /?^p^C;l] 

C2] /V^ffxC^plC 1 ; ] 

C3] SxZCCO. 5x/i/ + l] + (Z-^ZCiZ^,A:] )Cr 0 . 5X/I/ + 1] 

C4] RE*-RpO 
C5] Cfi’^CpO 
C 6 ] □•*‘5-<-+/ I , X 

C7] Z-(R,C)p(,X)U,X + i \R)o .xCpl^RG^i r/,I)’C/,-?] 

C 8] X-^Xr-iRI-^ + Z 0 . 5xZC;C0.5xr + l]+ZC;r0.5xC + l])o.xfpi 
C9] RE*-RE+RI 

CIO] Z^i)(r,/f )p( ,4?Z)C4 ,6}-? + (/?pl ) 0 . >< ( iC ) >‘ifG] 

Cll] A'-^’(ftpl) ".xCI^ + ZO. 5xZCC0.5xfl'+l; ] + ZCC0.5xi? + l;] 

Cl 2 ] CE-^CE+CI 

Cl3] -y6>^\EPSILON<l^ + / \ yXiS 

C14] Z*-X 

Cl5] -ylly\RESIDS*2 

C16] Z^6j(2,/V)p( ( ,/?6’o .xfff) iTV),,Z 

Cl7] ■*21^ iR0RMEFFECTS = 0 

Cl8] RE->-RE^RA-^(+/ RE) iR 

Cl9] CE*-CE^CA-(.^/CE)iC 

C20] TV->-TV^RA^CA 

C21] 'TYPICAL VALUE-. ' ;TV 

C22] 'ROty EFFECTS: ';RE 

C23] 'COLUMN EFFECTS: '-,CE 

V 
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VMSTATS^•J^SI 

V iV ;Z ;C iS irHit>jt:ANS iV AhS ;J 

L 1 1 r+-( is?z ) + . X 2 +-V- ( p^y' ) pi'jt'AI\/S-^{->-/-W ) + p ("2 1 1 , 1 , pV ) pVf pC^-^ ' ' 

[21 'BF8.2' Fw'I’( 0 , ipF) . [ 1 1 ( ipF) ,rf5o . x5^( [/y!l/^c’^ + /Z*2) * 

0.5 

[31 ^mEAt'iS iIj,6F8.2' FMT MEANS 

[41 'SSTD,DEV 3,BF8.2' FMT S^S i i N ^l) * 0 . 5 

[51 Z^\J^Q 

[6l ->-ENOx \{ qS) <J^J -^1 

[71 r^i[U^i^[ ;«/! 1 

[8l -*-6 ,Z+-Z,r'^0.5xr[[0.2 5x3i-,Vxx3l+F[[ 0,25^1i-'l/xx3l 

[9l ENDiZ^i ipS) ,3)pZ 

[10] 'iZlL. (^Fl”LFi!l,FF8.2 ' FMT Z[;ll 

[11] '^MEDIANS 0,FF8.2' FMT Z[;2] 

[12] 'DF. QFrL5ffl,FF8.2 ’ FMT Z[;3] 

V 



V/VFAiFFFV[D]V 

V NUMSUM X;N iQiTliCH 

[1] FF 

[2] ' NUMERICAL SUMMARY' 

[3] 70p'n» 

[4] iV^p;f^A'[U^(^^wiFF) / i^.a:] 

[5] 'SAMPLE SIZE = ' ; /!/ 

[6] 70p’r> 

[7] C-^0.5xA:[[0.125><7+A'xx7]fA[[ 0. 125x1 t-Fxi?] 

[8] 2'1^Q[4],0,F[4],0,0,( (Q[2l+F[6])i2),G[2].0 

[ 9 ] n^ri,«[6] , (e[6]nQ[2l),((F[l]^Q[7])T2) 

[101 ^^^.^[l] ,0,Q[71. (C[7lrQ[ll) , ( U[11+Z[F1 )t2) 

[111 T1*■T1,XL1^,0.XLN^,X^N'\■-X^1^ 

[ 12 ] Y’l-^ 4 5 pTl 

[131 CH^' MIDPTS , ,LOQ/ 8/MlN , JIEDIAN, ^UPQl8lMAX^ ^SPREADS , ' 
[14] 'FFl 0.2,0 I O' FMT Tl 
[151 70p’<t» 

V 



VFFFF[Q1V 
V H^ONEH XiAiBiC 
[ll /l^O.X.O 
[21 B*-X,0,0 

[31 r^O.O.Jf 
[4l 2+(r + fi'-r2x/3)J4 

[51 H*-Xll'\,2^H ,X\.QX^ 

7 
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vo^/ffaLG] V 

V Z^O/V&'3 A';yi;y2;J^l;Z2;/V;C 

[1] 71^(3xy[2])^2xA'[3] 

[2] 7 2^( 3xy[;V’l] )^2xy["2+,V^py] 

[3] 

[4] 72-^1 + 7,72 

[5] Z-^(yi,7,72)[( i/V) + l + (C = 7K7)’(7<72)=C-^71<72] 

V 



VP>iy?J’J/lL[G]V 

V P*-PARTIAL ^ \C \CIR\M\Y \X\BETA 
[Ij M^\1\qW 

[2] 'ENTER COLUMNS OP MATRIX TO BE PARTI ALLED OUT.' 

[3] /II :->/^2xiO=‘ + /'-(C-^G)€/Vj 

[4] Y*-Nl ;CIN‘>-(~MeC) /Ml 

[5] ^Q.pP-^CORRELATION R*-Y ^X + . ^ BET A-^Y\iX*-l , ( 2 + ( pV [ ; C ] ) , 1 ) [ ; C 3 

[6] A2:^Al,p'UNACCEPTABLE \/ A LUES, ENTER INTEGERS FROM 1 TO ' ‘,M 
. V 



V5P0VP£’5[G]V 

V Z*-SR0VRES X-,Y-,N\Q\D-,LW-,LF\U'M\UF 
Cl] N*pY-^YLiY*,Xl 

[2] ^-^0.5xy[L0.2 5x3 + .Vxi3]+7Cr 0.25xl + /l?xi3] 

[3] LN*-Qlll^D^QlZl^Q\.ll 

[4] UN->-QlZl+D 

[5] LF*-LN^Di2 

[6] UF^UN+Di2 

[7] y-(7>PF) + (7>Z/V) + (7>^;Cl] ) + (7>QC2] ) + (7>QC 3] ) + (7>i/V)+7>PP 

[8] Z-^'OOoooox;^ » [ 1+7] 

[9] Z^( (pZ[;l]),2xpZCl;])p(,Z) AND, ipZ)p' ' 

V 
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VRtiGRESSlQ^iV 

V Z-y REGRESS X. ■,R K vC -,K PXI fJV ',X PY BET A \ RSS \T SS -,S2 \E SS ^ I U j LE 

Cl] 2t( pz ) , 1 )p;r 

[2] X*-{Q,l-^EItrTERCEPT)\l,X 

C3] XPXlNV^'^i^X)* .y-X 

[4] BETA^XPXIMV+.^XPY^(.ISiX)*.^Y 

L5] RSS^i {Is^BETA)-^ .>^XPY)^C^{ ( + /Y ) * 2 ) wV^p , y 

[6] ESS-(TSS^( (iS)Y) + .xY)^c)-RSS 

C 7 ] S2^,ESSi iN-’D^K^ip , BETA )-METERCEPT 

[8] CR 

[9] ' AEOi/A' 

ClO] SOURCE, DF, SUM SQUARES, MEAN SQU A RE , F R AT 1 0 ' 

Cli] ’ ' 

[12] ’■l}REGRESSION\n,I^,BElB.R ' FMTiK ) , (.,RSS) , ( ,RSSiK) , ( ,RSSiK) iS^ 

[13] QR ^' ' 

[14] 'a fi’£5J/^6'^L.a,I4,tf£16. 4 ' FMT U U ^ 1) -K) , (, , ESS ) , S2 , 0 

[15] 'STOTAL a,J4,B£16.4' FMT i U ^ 1) , { ,T SS) , 0 ,0 

[16] ' ' 

[17] 'R SQUARE: •\,RSSrTSS 

[18] 'STD ERROR: ';,52*0.5 

[19] Cd*-'COEFFICIEUTS,T STATISTICS' 

[20] 'F15.4' f,'^7’i!(2,p,££T>^)p(,££T/3),(,££T/l)i(l 1 l^V^S2>^XFXI \\/ ) * 

0. 5 

[21] 'DO YOU UAUT A PRINTOUT OF TdE 1/ ARI ARC E ^C0\/ ART ARC E MATRIX?' 

[22] ->->11X1 'y ^IfQ 

[23] '\/ARIARCE^COVARIAUCE MATRIX: ',CH^'' 

[24] '£’12.4' FMT V 

[25] Al: ' DURBIR^Q ATSOU : ';(-^/((U,(T)-Cl + ,n)*2)v + /( ,C*-Y - X-^ . ^ bET A) * 

2 

[26] Z^^{2, U)pi., X^ .x-BETA) , ,C 

[27] B1:'D0 YOU UAUT TO FORECAST A VALUE FOR Y?' 

[ 28 ] ->-rixi »y ' 7=1+2 

[29] 'ENTER X VECTOR ('iK-,' VALUES)' 

[30] 'FORECAST OF Y VALUE: ' {C*- i El NT ERC EPT ) ^ 1 ,1.) ^ . x B ET A 

[31] 'VARIANCE OF FORECAST ERROR: ' -,S2xl^C ^ . xX P XI NV ^ . xlS)C 

[32] ->B1 

[33] C1:'D0 YOU N ART TO SCAT RESIDUALS VS. PREDICTED Y?' 

[34] ->-0x I ';i/ ' = 1 + 2 

[35] D£P^0. 5xyj/>-[/70, ( f / ( ( 0 . 7 5 x.V ) , 3 0 ) ) 

[36] SCAT Z 

V 
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75(7y5y’CmV 

V W^SCAT ZiNiXiYiCiHiUiSiLiliJ iKU'T;rLiGiDiBiAiC’,V 

[11 -►3x\(2 = + /2 = /^)v(x/W)> + /A/^pZ 

[2] Z^<S?(2,pZ)p( ipZ) .Z^.Z 

[3l y-*-Z [ ; 1 + iT-*- l + (pZ)[2l1 

[41 it^pZ^y^,Z[ ; 11 
[51 L-^i'’^5-^2pO 

[61 J^l-rQxp(D^NDIVX ,DlDIVY) ,B->-WID,DEP 
[71 ['2’'^10*[10®C’L^l£''‘2 0i-( U'LJ / Z) ^SLJ / Z) ^ DlJ 1 

[81 S^J^^UT^lS^J'\H'T^UT^li^CLfrUT^(l 2 5)x£/r1 
[91 ['[t7i-i'rxri'[./iTt'r 

[101 L[Jl-^l+Gx[ (B[J1«1) ^G^(i'[Jl«5[./1 ) H'T 

[111 z-^.r 

[121 + 1 
[131 /3^(<|>L)p0 

[I4l y-^l + [0.5 + (L[l1 -^l) X (y^5[11 ) fi'[ 11 1-5[11 
[151 y^l + [0.5 + (L[2l-l)x(rnB[2l)U'[2l"5[2l 
[161 I-^l 

[171 -►20xtl<r 

[181 /3[y[I;l1 ;y[I1Kl0L/3[y[7;l1 ;y[I11-i-l 

[191 -►18 + 6x/?<I-^I + l 

[201 J^l 

[211 O^0 = 7^;l[y[I;Jl ;I[I1 1 

[221 AlYLIiJl ;y[I1Xl0xl'>A:+l) + ( (A:+l)xA: = l/) + (A:^35t.2xc7)xB 
[231 -*-2lx + l 

[24l -*-21'x + 1 -*-1 

[251 t?-»-(<|)pi4)[iri-rL0.5-r(Lnl)x5^5"£' 

[261 /3[;C?[l11^/3[;C?[l11 + 36xO = /5[;C’[l11 

[271 /^[(?[2l;1-^/3[C[2l;1i-35xO = /3[^7[2l;1 

[281 023H56789miiKKJlIlBliGGFlEEDDCCBiAAf’\ ' Ll-^9Ai: 

1 

[291 'HANGE OF Xi ';5[l1,['[11 

[301 'RANGE OF Yt ';5[2l.t'[2l 

V 
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V5PLirCQ]V 

V Z^SPLIT X-,PiQ;RiSiT -,N il -,J iC 
Cn P-^'3iX 

[21 Q^~2il^X 
[3] R-<-~li2iX 
[41 5-^3+A' 

[51 I^( ( (P<Q)x/?>5)r(P>Q)xff<5)xQ = it 

[61 Q^iX^XLlI ,X ^XLNl) [ /I*-l ^\~3-<-R*-pX'\ 

[71 R^<-i2^Q)^2><X^I^ 

[81 p^;r[i^ 2 i 

[91 3x5'HA'[I-r4l )r-2xA'[I-*-5l 

[101 XII * 2'\^i R^1=J )■>■ iP^~ 1=J ) ->-Q>^0=J C = P<Q) n(Q<R)=C^P<R 
[111 ji:[z-t3l^(2’xi=.j)^(px"i=j)+pxo=j^(r=p<p) -(5<D=r-p<r 

[121 Z-^14"l+;f 

V 



VPry5I’IP7'IPP[aiV 

V STATISTICS RiXBAR-,VAR;l\l ;QiS ;Z 

[11 XBAR^(+/R)^!\I^q ,R 

[21 VAR-<-(-^/ URnXBAR)*2))iN^l 

[31 P^P[^P1 

[ 4 1 Q-^0 . 5 X c [ L 0 . 2 5 X 3 -r/l/ X 1 3 1 -rp [ r 0 . 2 5 X l+il/ X ; 3 ] 

[51 h^JEANx 'iXBAR 

[61 ^VARIAtJCEx 'iVAR 

[71 'STD, DEV.x ' \S^V AR*3.h 

[8l 'CCEFF. OF VARlATICNx ';Si\XBAR 

[91 'LOWER QL'ARTILEx ';Q[ll 

[10] 'UPPER QUARTILE: ';Q[3l 

[111 'MEDIAIU *;Q[2l 

[121 'TRIMEAIU ’ ;0.25xQ[ll-r^3[3>2xe[2l 
[131 'MIDMEARi ' ; ( ^ pZ ) x +/Z-^ ( ( P>S [ 1 1 ) a ( P <5 [ 3 1 ) ) /P 

[141 'RANGEx ';(r/P)-[/P 

[151 'MIDRAliGEx ' ; ( ( f /P ) ■<■[ /P) t 2 
[161 'MEAN ABSOLUTE DE.IATIONx ' ; ( + / 1 P cQ [ 2 1 ) ^ P 
[171 'INTERQUARTILE RANGE x ’;Q[3]-e[ll 

[181 'COEFE. OF SKEWNESSx ' ; ( { ^ { RnX B AR) * 3) i N 1) i S * 3 

[191 'COEFF, OF KURTOSISx ’ ; ( ( ( -r/ ( P-.A'P/IP ) * 4 ) i A/ - 1 ) r 5* u ) , 3 

7 
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1S'I:EMLEAF\.U^^ 

V Z^STEMLEAF A" ;T ; /? ;5 ; 51; I ; J ;F ; /5 ; ; L ; VV ; /J/3 ; 

[ 1 ] io®F^i5”2 0 + (A:[p;f ) cn ) i sc ale 

[2l 5I^(lf3x + /(Fxr)> 25 50)i- + /(pA’)> 25 100 

[31 Ar-L0.5-i-Jxrxl0* + /5I= 2 3 6 

[4] F->-^/ (SI=\9) / 0.5 2 1 1 0.5 2 2 1 0.5 

[51 /l-^-20p ' 0123456789/lBr55F55I.7 ' ' 

[6] I-^-FxLA[ n tIOxF 
[71 

[81 5^’n » [1+J^J[11>01 

[9 1 L-^l + ( |V-^-(A:<(10 xI)+Jx”i + 10xF)/A')7,10xL1I 
[101 AA^All^dO 10)t|[I1 
[111 ^{(pALL'\)^WW)/L1 

[121 Z-^Z,A'Wp5./^^, M (f/V + /l[Bl ).» + ', (Lir(p;J[Ll ) r,VV) ,2p ' ’ 

[131 -*-L2 

[141 LitZ-^z.A-;/p5,/5/i, ' I ' ,-4 [b1,a:;/p' ' 

[151 L2tI^I + fxi^(I = 0)x;f[ll<0 

[161 -*- 8 x\ 0 <pI-^( p {/) +1 

[171 Z^( ([ (pZ) tAT( 7) ,A’i7)pZ 

V 



75i'W5e[Dl V 

V SUMSQiC iNUMBERiBiBLDF 

[11 T^+^Y 

[21 i'’55'^(+/(+/Y*2) )tr^( {+/T)*2)iNUMBEH^+/R 

[31 BLSS^{->-/ ( {{B^->-/Y)*2) iK))r^C 

[41 MSBL^BLSSiBLDF^MAX’^1 
[51 ■*7^\{OPTICN = 2) 

[61 BLSS-^BLDF^Q 

[71 WTSS^TSS^BLSS*BTSS^(->-/ { (T*2)iR))-C 
[81 MSR^-BTSSiKr^l 
[91 TDF^NUMBERnl 

[lOl F^MSRiMSE^yTSSiEDF-^TDFr-BLDF^K^l 

[111 BLCCK^iBiK) nYBAR*-(+/T) ^NUMBER 

[121 » ANCVA TABLED 

[131 ' SOURCE DF SS MS f 

[141 ca-^" 

[151 'fl rFF>4!rW£A/?’(!).I5.F13.2,F11.2,F8.2' FMT (K^.!) ,BTSS ,MSR , F 
[161 ^(CP’riON = l) /L5 

[171 'Q BLFFF51!1,I8,F13.2,F11.2,F8.2 ' FMT ( M AX -1) , B LS S ,MS B L , {MSB LiMS £ 

) 

[181 L5t’Q FFF5F1!1,I9,F13.2,F11.2' FMT EOF ,WTSS J^SE 

[191 ’ID rt)7’/5LI!l,l9,F13.2> FMT TDF^TSS 

[201 'm’^SQUARE = ffl,F5.3» FMT {BTSS-^BLSS) iTSS 

[211 '\!}OVERALL MEAN = D.F10.2' FMT YBAR 

[221 'aTREATMENT EFFECTS 1I1.F6.2' FMT {T i N) -^Y B AR 

[231 -*{0PTI0N = 1) /O 

[241 '71BLCCK EFFECTS ffl,F6.2' FMT BLOCK 

V 
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vry/CttCDlv 

V A^TiVICt. x-,y 

[1] Z-Y + A3RSR X-^Y*-A3RSR X 

V 



^UTCORDLBly 

V UTCOND X 

[1] X-^i2pipX) ,1) pX 

[2] Z-'» 

L3] I^Ox/l/ff^(p;f )[2] 

[4] i^-py-(7C^y^(y^MI56’)/y^A [ ;7+l] ] 

[5] (^-^0. 5xy[L0. 125x7+.Vxi7]+i/Lr 0 . 1 2 5 x 1 +,V x i 7 ] 

[6] -<-7+2xLA'y;’.n2 

[7] Z^Z,2C4] 

[8] ->4 + l lx/V7? = 1 

L9] Z-Z,5C4], i^[6]-(i[2] 

[ 10 ] -^-4 x;i/it> 7 ^ 7+ 1 

[11] Z-Z,y[l], Q[2],i^[4], Q[6] ,y[/'y], V 

[12] ^4 + 1 3x/;//f = j^J+ 1 

[13] Z^Z,y[l]w[l 2 4 6 ll,UlR'\J\l 

[14] ^4 + 16x/'/ii = J^J+l 

[15] CR^' MEDIANS' 

[16] ^0xpa^»F10.2 ' FhT Z^(RR,l)pZ 

[17] Z-(.'i;/^, 6)pZ 

[18] CH *- ' MIN , LOQ , MEDIAN,UPw , -lAX , SIZE ' 

[19] ^-OxpiJ-'f'lO. 2 ' FMT Z 

[20] Z-^(,V/i , 8) pZ 

[21] Ca-^'MIN ,LOS, LOQ , MEDIAN ,UPQ ,UPS, nAX , SIZE ' 

[22] U^'^’9.2' FMT Z 

V 
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